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I. INTRODUCTION 



In the next year, the Large Hadron CoUider (LHC) will begin operation at CERN, usher- 
ing in a new regime of directly probing physics at the shortest distance scales. The LHC will 
search for physics beyond the Standard Model, as well as for the Higgs boson. The Higgs 
mechanism, which describes the breaking of electroweak symmetry in the Standard Model 
and its supersymmetric extensions, is simultaneously the keystone of the Standard Model 
(SM) and its least-well-tested ingredient. 

The dominant process for Higgs-boson production, over the entire range of Higgs masses 
relevant for the LHC, is the gluon fusion process, gg H, which is mediated by a heavy- 
quark loop [1]. The leading contribution comes from the top quark. The contributions from 
other quarks are suppressed by at least a factor of (9(m^/mj), where rrit, rrib are the mass of 
the top and of the bottom quark, respectively. Because the Higgs boson is produced via a 
heavy-quark loop, the calculation of the production rate is quite involved, even at leading 
order in a^. The inclusive production rate ior gg ^ H + X has been computed at next-to- 
leading order (NLO) in [2], including the full quark-mass dependence [3], which required 
an evaluation at two-loop accuracy. The NLO QCD corrections increase the production rate 
by close to 100%. However, in the large-rrit limit, namely if the Higgs mass mn is smaller 
than the threshold for the creation of a top-quark pair, ttih < 2mt, the coupling of the 
Higgs to the gluons via a top-quark loop can be replaced by an effective coupling [2, 4, 5]. 
This approximation simplifies calculations tremendously, because it effectively reduces the 
number of loops in a given diagram by one. 

It has been shown that one can approximate the full NLO QCD corrections quite accu- 
rately by computing the NLO QCD correction factor, K^^*^ = a^^'^/a^'^, in the large-m^ 
limit, and multiplying it by the exact leading-order calculation. This approximation is good 
to within 10% in the entire Higgs-mass range at the LHC, i.e. up to 1 TeV [6]. The reason 
it works so well is that the QCD corrections to gg ^ H are dominated by soft-gluon effects, 
which do not resolve the top-quark loop that mediates the coupling of the Higgs boson to 
the gluons. The next-to-next-to-leading order (NNLO) corrections to the production rate 
for gg ^ H have been evaluated in the large-m^ limit [7] and display a modest increase, less 
than 20%, with respect to the NLO evaluation. The dominant part of the NNLO corrections 
comes from virtual, soft, and coUinear gluon radiation [8], in agreement with the observa- 
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tions at NLO. In addition, the threshold resummation of soft-gluon effects [6, 9] enhances the 
NNLO result by less than 10%, showing that the calculation has largely stabilized by NNLO. 
The soft and coUinear terms have recently been evaluated to one more order (N^LO) [10], 
with partial results to N^LO [11], which further reduces the uncertainty on the inclusive 
production cross section. 

Backgrounds from SM physics are usually quite large and hamper the Higgs-boson search. 
A process that promises to have a more amenable background is Higgs production in as- 
sociation with a high transverse-energy (Et) jet, pp ^ H + jet. In addition, this process 
offers the advantage of being more flexible in the choice of acceptance cuts to curb the back- 
ground. The pp ^ H + jet process is known exactly at leading order [12], while the NLO 
calculation [13] has been performed in the large-m^ limit. For Higgs -|- jet production, the 
large-mj limit is valid as long as uih < 2 and the transverse energy Et is smaller than 
the top-quark mass, Et < rrit [14]. As long as these conditions hold, the jet-Higgs invariant 
mass can be taken larger than rrit [15]. At itlh = 120 GeV, the NLO corrections to the 
pp ^ H + jet process increase the leading-order prediction by about 60%, and are thus of 
the same order as the NLO corrections to fully inclusive pp — > if production considered 
above. At present, the NNLO corrections to pp ^ H + jet are not known. 

An even more interesting process is Higgs production in association with two jets. In 
fact, a key component of the program to measure the Higgs-boson couplings at the LHC is 
the vector-boson fusion (VBF) process, qq — > qqH via t-channel W ot Z exchange. This 
process is characterized by two forward quark jets [16]. The NLO corrections to Higgs 
production via VBF fusion in association with two jets are known to be small [17]. Thus, 
small theoretical uncertainties are predicted for this production process. The production of 
H + 2 jets via gluon fusion is a part of the inclusive Higgs production signal. However, it 
constitutes a background when trying to isolate the HWW and HZZ couplings responsible 
for the VBF process. A precise description of this background is needed in order to separate 
the two major sources oi H + 2 jets events: one needs to find characteristic distributions 
that distinguish VBF from gluon fusion. 

The production oi H -\-2 jets via gluon fusion is known at leading order in the large-mj 
hmit [18] and exactly [19]. As in the case ol H -\- 1 jet production, the large-m* hmit is 
valid as long as ttih < 2 rrit and the jet transverse energies are smaller than the top-quark 
mass, Et^i,Et^2 < ^t, even if the dijet mass, as well as either of the jet-Higgs masses, is 
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larger than the top-quark mass [15]. At present, the NLO calculation for the production 
of if + 2 jets via gluon fusion has not yet been completed/ not even in the large-mt limit, 
although the necessary amplitudes, namely the one-loop amplitudes for a Higgs boson plus 
four partons [21] and the tree amplitudes for a Higgs boson plus five partons [22, 23, 24], 
have been computed. Considering the aforementioned results for the NLO corrections to 
fully inclusive Higgs production and to Higgs production in association with one jet, and the 
fact that the largest NLO corrections are usually found in gluon-initiated processes, there 
is no reason to expect that the NLO corrections to H + 2 jets production via gluon fusion 
are small. Thus, such a NLO calculation is highly desirable. 

Another distinguishing feature of VBF is that to leading order, and with a good approxi- 
mation also to NLO, no color is exchanged in the t channel [25]. The different gluon-radiation 
pattern expected for Higgs production via VBF versus its major backgrounds, namely tt pro- 
duction and QCD WW + 2 jet production, is at the core of the central-jet veto proposal, 
both for heavy [26] and light [27] Higgs searches. A veto of any additional jet activity in the 
central-rapidity region is expected to suppress the backgrounds more than the signal, be- 
cause the QCD backgrounds are characterized by quark or gluon exchange in the t channel. 
The exchange of colored partons should lead to more central gluon radiation. In the case of 
Higgs production via gluon fusion, with two jets separated by a large rapidity interval, the 
scattering process is dominated by gluon exchange in the t channel. Thus, as for the QCD 
backgrounds, the brcmsstrahlung radiation is expected to occur everywhere in rapidity, and 
Higgs production via gluon fusion can be likewise checked by requiring a central-jet veto. 
For hard radiation, the effectiveness of such a veto may be analyzed through H + 3 jets 
production, which for gluon fusion is known in the large-mt limit at leading order [22]. An 
analogous study that also includes soft radiation has just been completed, by interfacing 
matrix-element calculations oi H + 2 and 3 jets production with parton-shower effects [28]. 

The discussion above motivates us to look for methods to compute Higgs production 
in association with many jets at least to NLO or higher accuracy. Tree amplitudes for a 
Higgs boson produced via gluon fusion together with many partons have been computed 
analytically in the large-m^ limit [23, 24] and are available numerically through the matrix- 
element Monte Carlo generators ALPGEN [29] and MADEVENT [30]. However, at one 

^ Preliminary NLO results have been reported recently [20] . 
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loop only the QCD amplitudes for a Higgs boson plus up to four partons are known [21]. 
The three-parton helicity amplitudes were computed analytically in the large- limit [31]. 
The four-parton results are semi- numerical, except for two cases: (1) analytic results have 
been presented for the case of two quark- antiquark pairs [21], and (2) the case of four gluons 
all with the same helicity is now known analytically as well [32] . 

Recently a lot of progress has been made in the computation of tree-level gauge-theory 
amplitudes [23, 24, 33, 34, 35, 36, 37, 38, 39, 40, 41], following Witten's proposal [42] 
of a weak-weak coupling duality between = 4 supersymmetric gauge theory and the 
topological open-string B model in twistor space (for a review, see ref. [43]). Recursion 
relations for computing tree amplitudes have been written [36, 37] which employ only on-shell 
amplitudes at complex values of the external momenta. Likewise, progress has been made in 
the computation of one-loop amplitudes [44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57], 
particularly those with many external legs, for which there are many different kinematic 
variables. If careful attention is not paid to the analytic structure, very large analytic 
expressions may result. One approach to this problem is to compute the loop amplitudes 
numerically or semi-numerically [20, 21, 58]. Another approach is to exploit the analytic 
properties of the amplitudes in order to facilitate their computation. The two main analytic 
properties of loop amplitudes are branch cuts and factorization poles. 

In pure QCD with massless quarks, there are special helicity configurations for which 
the tree amplitudes vanish, and the corresponding one-loop amplitudes are finite. Such he- 
licity amplitudes are of no immediate use in phenomenology, because they contribute first 
at NNLO. However, they have interesting analytic properties. The specific helicity ampli- 
tudes in this category are the pure-glue one-loop amplitudes for which all the gluons have 
the same helicity, or else all but one do, namely, An\i^, 2"*", 3+, . . . , and in addition 
the amplitudes with one pair of massless external quarks and {n — 2) positive-helicity glu- 
ons, ""(l^, 2+, 3^, . . . , n"*"). These amplitudes contain no branch cuts; they are purely 
rational functions of the kinematic variables [59, 60]. Compact forms for the amplitudes 
have been found by constructing on-shell recursion relations, along the lines of the tree-level 
relations [36, 37], and then solving them in closed form [52, 53]. 

All the other one-loop hehcity amplitudes in QCD contain branch cuts. The terms with 
branch cuts can be determined from unitarity in four dimensions; then recursion relations can 
be established for the remaining rational terms [54, 56] . This approach offers an efficient way 
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in principle of building up a generic n-point one-loop amplitude from one-loop amplitudes 
with a smaller number of points. It has been applied in practice to determine the rational 
parts of the n-gluon amplitudes with two negative-helicity gluons in arbitrary locations in the 
color ordering [55, 57], and also for cases with three or four color- adjacent negative-helicity 
gluons [56]. (For n — 6, the rational terms for three color- nonadjacent negative-hehcity 
gluons, as well as two negative-helicity ones, have been worked out [61] using Feynman- 
diagrammatic methods [62].) 

A closely-related framework can be constructed for amplitudes for a scalar field, such as 
the Higgs boson H, interacting with an arbitrary number of quarks and gluons. In the large- 
rrit hmit, the Higgs field couples to gluons via the effective operator H tr G jj,i,G'^'^ induced 
by the heavy-quark loop. To make the situation parallel to the pure-QCD case, we follow 
refs. [23, 24] and treat H as the real part of a complex field cf). Here couples to the self- 
dual part of the gluon field strength via the operator 0tr Gsd fiuGso, whereas the conjugate 
field (f)'^ couples to the anti-self-dual part. Then, exactly as in the case of pure QCD [i.e. 
when is absent), the amplitudes An{(p, 1^, 2+, 3"*", . . . , n"*") and An{4>, 1^, 2+, 3"*", . . . , n"*") 
vanish at tree level, and are free of branch cuts at one loop. (They are also free of infrared 
and ultraviolet divergences.) Thus the one-loop amplitudes, An\(p, 1^, 2+, 3+, . . . , n+) and 
An\(f), 1^, 2^, 3"*", . . . , 71+) , can be written as rational functions of the kinematic variables. 
The 0^ amplitudes are obtained from the amplitudes by parity, which also reverses all 
parton helicities. Amplitudes containing a scalar Higgs field H are given by the sum of 
the 4> and 4>^ amplitudes; whereas amplitudes containing a pseudoscalar field A instead are 
found from the difference of the two amplitudes. In these sums and differences, the finite 
(f) amplitudes have to be combined with the image under parity of divergent amplitudes 
with mostly negative gluon helicities. 

The aim of this paper is to construct and solve on-shell recursion relations for the finite 
one-loop amplitudes with a single scalar field 0. To assemble complete analytic helicity am- 
plitudes for H or A fields, other, divergent helicity configurations are also required. These are 
the (f)^ amplitudes with the same helicity configurations as those computed here, or equiva- 
lently, the amphtudes with all parton helicities reversed. The full set of hehcity amplitudes 
for H or A fields requires as well amplitudes with two or more negative parton helicities 
and two or more positive helicities. We leave such computations to future work — except 
that we shall quote the recent results of ref. [32] for four identical-helicity gluons. We stress. 
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however, that the amphtudes we compute represent one piece in the (/>-(/)''' decomposition of 
H or A amplitudes for which the corresponding tree amphtudes arc nonvanishing. Hence 
they contribute to NLO Higgs cross sections. (This is in contrast to the finite amphtudes 
in pure QCD, without a field, which first contribute at NNLO.) A second reason why the 
finite (p amplitudes are useful is that in the recursive approach we pursue, amplitudes with 
few negative helicities are needed to supply factorization information for amplitudes with 
more legs and more negative helicities. 

The construction of on-shell recursion relations relies on information about the factor- 
ization properties of amplitudes, as momentum invariants built from two or more external 
momenta become null. In the "coUinear" case of two external momenta, this behavior, in- 
volving complex momenta, is quite subtle at the loop level. Loop amplitudes can contain 
double poles in coUinear invariants [52], and "unreal" poles for which the coUinear limit in 
real momenta is not singular, yet nonetheless a single pole arises for complex momenta [53] . 
In our recursive approach to the finite quark-gluon amplitudes, An\(j), 1^, 2+, 3+, . . . , 
we encounter both these types of poles. A knowledge of the residues at these poles is essential 
for constructing correct recursion relations. It might be possible to rigorously analyze such 
contributions using "space-cone" gauge techniques [63] . In this paper we adopt a much more 
pragmatic approach; indeed, we argue that the unreal poles take precisely the same form as 
in the pure-QCD case [53], for which several types of consistency checks were available. 

Using this approach, we present recursion relations for the finite 0-quark amplitudes with 
an arbitrary number of external gluons. We find compact solutions to the relations, valid for 
all n. To confirm these relations we perform nontrivial consistency checks of the factorization 
properties of the solutions. One of the factorizations, as the momenta of the quark and anti- 
quark become coUinear, is onto the one-loop 0-plus-n-gluon amplitudes with a single negative 
gluon helicity and the rest positive, and determines this sequence of amplitudes. 

The 0-amplitudes with two or more negative helicities can be obtained by extensions of 
the methods of refs. [47, 54, 56, 57, 64, 65]. Unitarity (or the related loop-level application of 
maximally- helicity- violating (MHV) rules [44]) can be used to determine the cut-containing 
functions from known tree- level amphtudes [47, 64, 65]. The cut-containing parts of Higgs 
amplitudes with an arbitrary number of identical-helicity gluons have been found in this 
way [32]. The remaining rational functions can then be computed via a factorization boot- 
strap approach, in the form of on-shell recursion relations analogous to the approach of 
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ref. [54, 56, 57] in the pure-QCD case. In this pubhcation, however, we mainly concentrate 
on the amphtudes with one or no negative hehcitics, which do not contain any cuts. 

The paper is organized as follows: After introducing some convenient notation in sec- 
tion II, we review the on-shell recursion relations in section III and the previously-known 
one-loop finite amplitudes in pure QCD in section IV A. In section IV B we summarize known 
results for amplitudes with a cf) field and present a recursion relation for the MHV tree-level 
0-plus-n-gluon amplitudes, An\4>, l"*", . . . , i~, . . . , j~, . . . , n+). In section V we discuss some 
interesting features of the soft-Higgs limit of the 4> amplitudes, as the momentum of van- 
ishes. In the tree-level case, the soft-Higgs limit always produces a fixed, helicity-dependent 
multiple of the corresponding QCD amplitude. At one loop, however, this is not true in 
general, although it is true for the finite (f) amphtudes [23] . 

In section VI we present and solve the recursion relation for the one-loop 0-plus-n-gluon 
amphtude, with all the gluons of like helicity. In section VII we do the same for the one- 
loop 0-quark-gluon amplitude An\(f), 1^, 2+, 3^, . . . , n"*"). Using factorization, we obtain the 
0-plus-n-gluon amplitudes with a single negative-helicity gluon, An \(f), 1^, 2+, 3+, . . . , n"*"). 
Section VIII collects explicit results for all helicity configurations for the 0-1-2 and 3 parton 
amplitudes, as well as partial results for 0-1-4 partons, using also results from refs. [31, 32]. 
In section IX we draw our conclusions. In appendix A we discuss the normalization of the 
one-loop amplitudes. In appendix B we prove the validity of the on-shell recursion relations 
for 0-amplitudes at tree level. 

II. NOTATION 

In this section we first discuss the decomposition of amplitudes containing a single Higgs 
boson and multiple partons, into amplitudes with a and a 0^ field, respectively. Then we 
review the color-decomposition of these amplitudes at tree level and at one loop. Finally, 
we introduce some convenient notation for the manipulation of these amplitudes in the 
spinor-helicity formalism. 
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A. (j)-(t>^ decomposition of Higgs amplitudes 



In ref. [23], the MHV rules for pure-gluon amplitudes [33] were extended to include the 
coupling to the SM Higgs boson in the large-mf limit. These rules were then further extended 
to compute amplitudes with quarks [24] . Although we pursue in this paper a recursive rather 
than MHV approach, the construction of refs. [23, 24] is still useful, so we briefly summarize 
it here. 

The couphng of the SM Higgs boson to gluons is through a fermion loop [1, 4, 66], with 
the dominant contribution coming from the top quark, because the Higgs coupling to quarks 
is proportional to the respective quark masses. For large top mass, nit oo, the top quark 
can be integrated out, yielding the following effective interaction [4, 5] 

£'^* = ^//trG'^.G^^ (2.1) 

where the strength of the interaction is given, to leading order in the strong coupling, by 
C = as/{6TTv), with v = 246 GeV. (The trace is normalized so that trC^^G"^ = G^^G"'''^.) 

In ref. [23], it was found that the MHV or twistor-space structure is simplest upon 
dividing the effective interaction Lagrangian into two parts, a holomorphic (self-dual) and 
an antiholomorphic (anti-self-dual) part. Then the MHV rules for QCD [33] can be extended 
straightforwardly to include these new interaction types. The Higgs boson is considered to 
be the real part of a complex field, = j{H + iA), so that 

C 



H,A - Y 

= C 



>tv G SD ixuGgjy -\- (f)^ Gasdhi^G'^qj-, 



(2.2) 
(2.3) 



Here the field strength has been divided into a self-dual (SD) and an anti-self-dual (ASD) 
field strength, 

GZ = liG^"" + , GZo = ^(C'''' - *G'"') , *G^"' = y'"'f"'Gp^ . (2.4) 

The scalar H and pseudoscalar A are reconstructed from the complex fields (f) and 0^ ac- 
cording to, 

H = (P + (P\ A = -{(P-(P^). (2.5) 

Prom eq. (2.5), the amplitude for a single scalar Higgs boson plus multiple quarks and gluons 
can be recovered, at any loop order I, as the sum of the amplitudes with and 0^, 

A^\H,...) = Al^\<P,...)+Al^\<P\...), (2.6) 
9 



where ". . ." indicates an arbitrary configuration of partons. 

As a byproduct, it is trivial to obtain tlic amplitudes for a pscudoscalar A plus partons 
as well, in the limit where the coupling is described by the effective Lagrangian (2.2), 



When using eq. (2.7) one should keep in mind that the value of the normalization factor C 
in the pseudoscalar case differs from that in the scalar case. For example, if the pscudoscalar 
state arises from a two-Higgs doublet model, and the only surviving contribution is from the 
top quark, with rrit taken to be large, then the leading-order value is C — ag cot /3 / {Attv) , 
where tan/3 = V2/V1 is the ratio of Higgs vacuum expectation values. In contrast to the 
scalar case, this coupling does not get renormalized by QCD corrections [67, 68], although 
at order the pseudoscalar A begins to couple to the divergence of the light-quark axial 
current [69]. As we shall see in section IV C, the amplitudes produced by this operator 
vanish in the limit that the masses of the light quarks go to zero. So there is no contribution 
to the cross section at NLO, that is, at order C as ~ a^. 

The amplitudes with are related to those with 0^ by parity. 



That is, to go from an amplitude with to an amplitude with (f)^ one needs to reverse the 
helicities of all gluons, and replace {ij) with [j i]. It is therefore sufficient to compute only 
the (f) amplitudes, and get the (f)^ amplitudes by parity. Note that reconstruction of the 
scalar H and pseudoscalar A helicity amplitudes (2.6) and (2.7) from amplitudes, with 
the aid of eq. (2.8), requires pairs of (f) amplitudes with reversed parton helicities. 

B. Tree- level color decompositions 

We are interested in calculating tree-level and one-loop amplitudes for a single color- 
neutral scalar field, 0, plus either n gluons or two quarks and {n — 2)-gluons. Because 
is color neutral, the color organization of these amplitudes is identical to that of the 
corresponding amplitudes in pure QCD. We now review the standard, trace-based color 
decompositions [70, 71, 72, 73, 74] for QCD, but add the argument to all the amplitudes. 
By "QCD", we actually mean a slight generalization: SU(A^c) gauge theory with Uf massless 



A^\A,...)=^^[A^\<P,...)-A^\<P\...)] . 



(2.7) 




(2.8) 
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quarks (fermions in the fundamental Nc + Nc representation) and Ug massless squarks 
(complex scalars in the Nc + Nc representation). 

In general, the coefficients of the various color-trace structures, called partial amplitudes, 
are built out of several primitive amplitudes [74]. Primitive amplitudes are color-ordered 
building blocks, which are functions only of the kinematic variables. We will write recursion 
relations for the different primitive amplitudes. 

For tree amplitudes containing the field (j) and n external gluons, the color decomposition 
is [23, 70, 71, 72], 

= Cg'^-'^ 2^ Tr(TMi)TM2)...TMn)) (0^(^(1/11^ 2''^...,r^''")), (2.9) 

where Sn is the full permutation group on n elements, Zn is the cyclic subgroup preserving 
the trace, and j^^ denotes the j-th (outgoing) momentum kj and helicity hj. The are 
fundamental representation SU(A^c) color matrices normalized so that Tr(T"T'') = 5"*. We 
only need to compute the partial amplitudes >ln '*(0, 1^*% 2^*2^ ^ ^hn-^ having the standard 
cyclic ordering; the remaining quantities entering eq. (2.9) are obtained by applying the 
permutations a to the momentum labels. (Whenever the permutation a acts on a list of 
indices, it is understood to be applied to each index separately: cr(3, . . . , n) = o-(3), . . . , (j(n), 
etc.) The normalization factor C is to be set to ag/iQiiv) for the SM Higgs boson in the 
large-mt limit, and to 1 in the pure- Q CD case. 

For tree amplitudes containing a (j) field, a quark pair in the fundamental representation, 
and {n — 2) external gluons, the color decomposition is [72] , 

^(°)(0,l„2„3,...,n) = Cg^-^ (T'^^^^' • • • ^"''^"04^ 4°H</', 1,-, 2,; a(3, . . . , n)) . 

(2.10) 

Here we have suppressed the helicity labels, and Sn-2 is the permutation group associated 
with the n — 2 gluons. Because the color indices have been removed, there is no need 
to distinguish a quark leg q from an anti-quark leg q in the partial amplitudes; charge 
conjugation relates the two choices. Helicity conservation implies that the helicities of the 
fermionic legs 1 and 2 are opposite. We take the hehcity of leg 1 to be negative, and that of 
leg 2 to be positive. The other case is obtained by using a reflection symmetry, 

4°)(0, 1+, 27, 3, 4, . . . , n) = (-1)MW(0, 27, 1+, n, n - 1, . . . , 3) , (2.11) 
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and then relabelling the external legs. 



C. One- loop color decompositions 

At one loop, the color decomposition for a (p field plus n external gluons contains double 
traces as well as single traces [73] , 

K2J+1 

^W((/),l,2,...,n) = C^"cr Yl Yl Gr„;» ^^(l, 2, . . . , n)) , (2.12) 

where \_x\ is the largest integer less than or equal to x. The leading-color structure, 

Gr„;i(l) = iVe Tr(T"i---T'^"), (2.13) 

is Nc times the tree color structure. The subleading-color structures are given by 

Gr„;c(l) = Tr(r"i • • -T"^-!) Tr(r"^ • • -T"") . (2.14) 

In eq. (2.12), Sn-c is the subgroup of Sn that leaves Gr„;c invariant. We have extracted a 
loop factor, relative to ref. [73], of 

1 r(l + e)r2(l-e) 1 , 

For the finite helicity amplitudes, which form the main subject of this paper, cp may be set 
to l/(47r)^. The A^-i are the "leading-color" partial amplitudes, while the A^-c for c > 1 
are subleading color, because for large An-i alone gives the leading contribution to the 
color-summed correction to the cross section, obtained by interfering An^ with An^ . The 
0{as) corrections to the operator coefficient C [75] (see appendix A) can also be included 
in eq. (2.12). 

For general helicity configurations, A^.^i has the following dependence on the numbers of 
fermions and scalars, Uf and n^, 

where ^If''^'*^ gives the contribution of a (gluon, fermion, scalar) in the loop [73]. For the 
finite helicity configurations, supcrsymmctry Ward identities [76] — valid for the pure-QCD 
case only — imply that a!}^ — —^v) = -^n-, or in other words that, 

A.,i(l±, 2+, . . . , n+) = (l - ^ + ^)4^1(l^ 2+, . . . , n+) . (2.17) 

12 



For the finite helicity amplitudes containing a single (j) field and n gluons, we shall see in 
section VI that the fermion and scalar loop still have equal and opposite contributions, 
An"^ — —J'^n- However, the gluonic loop differs, so the equation analogous to eq. (2.17) is, 

A„;i(0,l±,2+,...,n+) = A|f-^](0,l±,2+,...,n+) 

where 

= 4^1 - . (2.19) 

We emphasize that the relation A[f^(0, . . .) = — An'(0, . . .), i.e. the decomposition used in 
eq. (2.18), is only valid for the finite helicity amplitudes. 

The fermion and scalar loop contributions, proportional to n/ and respectively, only 
enter the leading-color partial amplitudes A^-x- The subleading-color contributions to the 
one- loop n-gluon amplitudes, A„;c for c > 2, come just from purely-gluonic graphs. This 
result holds in pure QCD, and for amplitudes containing a field. The subleading-color 
terms are given by sums over permutations of the gluonic contributions to the leading-color 
terms [64]. The formula is, 

^„;e(0,l,2,...,c-l;c,c+l,...,n) = (-1)'=-^ J] A^^\^,a{\,2, . . . ,n)) , (2.20) 

(Tecop{Q}{/3} 

where G {a} = {c- 1, c-2, . . . , 2, 1}, /5i G {/?} = {c, c+ 1, . . . , n-1, n}, and COP{a\{fi} 
is the set of all permutations of {1,2, with n held fixed that preserve the cyclic 

ordering of the ctj within {a}, and of the /3i within while allowing for all possible 
relative orderings of the ctj with respect to the /Sj. Equation (2.20) was established in the 
pure-QCD case, but the arguments rely just on the property that the color structures of 
the one-loop graphs involve only structure constants . This property still holds when 
the vertices from the interaction (p ti G s d hi^G's'd ai'c included. For another version of the 
argument, which uses /"'"^-based color structures instead of trace-based ones, and in addition 
avoids the introduction of the subleading-color pieces A^-c: see ref. [77]. 

The color decomposition at one loop, for amplitudes containing a (p field, a pair of fun- 
damental representation quarks, and {n — 2) gluons, is [74] , 

n-1 

^W(0,l„2„3,...,n) = Cg-crJ2 E Grg)(<^(3, ■ ■ ■ , n)) A„;,-(0, 1„ 2,; <j(3, . . . , n)) , 

j=l aeS„-2/Sn;j 

(2.21) 
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(b) 



FIG. 1: (a) An L type primitive amplitude, in which the fermion hne turns left on entering the 
loop (following the arrow), (b) In an R type primitive amplitude, the fermion line turns right. 





FIG. 2: (a) Graphs in which the external fermion line passes to the left of the loop are assigned to L 
type, (b) Graphs in which it passes to the right are called R type. Gluons, fermions or scalars run 
in the loop. The same decomposition is used when gluons are emitted from the external fermion 
line. 

where we have again extracted an extra factor of cr, and the color structures Gr^^^^ are 
defined by. 



Grgf(3,. 
Grg(3;4,. 

,j + i;j + 2,. 

'^^n;n-ll'-'5 • 



112 ' 



.,n) 
. ,n) 



0, 



T--U\ j = 3,...,n-2, 



'^2 



Tr(T 



(2.22) 



As in the n-gluon case, the An-^i are leading-color partial amplitudes, and the An-^j for j > 1 
are subleading color. 

Following ref. [74], where more details can be found, we introduce the primitive ampli- 
tudes, 

A^^^lf(/),lf,3,4,...,?,2,, ? + l,...,n), 

(2.23) 



^P](0,V,3,4,...,j,2^,j + l, 



J = 1,|,0, 



where J = \ and J = denote the contributions with a closed fermion loop and closed 
complex scalar loop, respectively, and "L" and "i?" denote the two orientations of the 
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fermion line in the loop. Generic diagrams contributing to the L and R terms are shown 
in figs. 1 and 2. Because the primitive amplitudes can be used to build amplitudes with 
any color representation for the fermions, we label them by / to denote a generic fermion in 
any color representation. Diagrams without closed fermion or scalar loops are assigned to 
J = 1; they may or may not contain a closed gluon loop, as the two types of diagrams mix 
under gauge transformations. For notational simplicity, we suppress the superscript "[1]", 
= An''^^ , = An'^^\ The primitive amplitudes (2.23) are not all independent. The set 
of diagrams where the incoming leg 1 turns left is related, up to a sign, to a corresponding 
set where it turns right. Thus, the two sets are related by a refiection which reverses the 
cyclic ordering, 

V, 3, 4, . . . , 2y, . . . , n - 1, n) = ly, n, n - 1, . . . , 2;, . . . , 4, 3). (2.24) 

The leading-color contribution to eq. (2.21), An-i, is given in terms of primitive amplitudes 

by, 

A„.i(0, 1^, 2,; 3, . . . , n) = A^icf^, 1/, 2^, 3, . . . , n) - ^A^(0, 1/, 2^, 3, . . . , n) 

c 

+ ^>l^'[^/^l(0,l/,2/,3,...,rz) + ^>l^'[°l(0,V,2;,3,...,rz). (2.25) 

For QCD the number of scalars vanishes, Hs = 0, while nf is the number of light quark flavors. 
The subleading-color partial amplitudes A„y>i appearing in eq. (2.21) may be expressed as 
a permutation sum over primitive amplitudes, 

An;j{(l>, Ig, 2g; 3, . . . ,i + 1; J + 2, j + 3, . . . , n) 



<TeCOP{a}{t3} >- 



A^'W(0,^(l/,2/,3,...,n)) 



^^y^ 2,, 3, ... , n)) - ^^^'M (0, a{y, 2,, 3, ... , n))J , (2.26) 

where e {a} = {j + 1, j, . . . , 4, 3}, /3j e {f3} = {1,2, j + 2,j + 3,...,n - l,n}, and 
COP{a}{P} is the set of all permutations of {1,2, . . . ,n} defined after cq. (2.20), except 
that here leg 1 is held fixed. However, as in the n-gluon case above, using /"^'^-based color 
structures [77] instead of trace-based ones, it is possible to avoid the introduction of the 
subleading-color amplitudes An-j>i. 

In the case studied below, where all external gluons carry the same helicity, the fermion 
loop and scalar loop are the same up to a sign, 

^^'[°1(0, 17, 2\...,j;,..., n+) = -^^'1^/^1(0, 17, 2+, . . . , j/, . . . , n+) 
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= A^(0,l7,2+,...,j/,...,n+). (2.27) 

As mentioned above, in contrast to pure QCD, this relation does not follow from supersym- 
metry; instead we will show it recursively below. Then by computing the closed scalar-loop 
primitive amplitude we obtain also the closed fermion-loop primitive amplitude. 
We shall find it convenient to compute the combinations 

and A^^^A^-A^ (2.28) 

instead of and A^. Also in contrast to pure QCD, for which an identity re- 
lates A^^*(ly^, 2+, . . . , jj^, . . . , n"*") with the amplitude with the cyclic ordering reversed, 
A^~^{1J , . . . . . . , 2"*") [53], here we will have to compute all for j from 2 up 

to n. The generic recursion relation will only hold for j < n — 2. However, we will be able 
to compute the cases of A^'^ ior j — n or j — n—1 by taking coUinear limits of appropriate 
amplitudes with larger n. The scalar-loop contribution A^ vanishes iij — norj — n— 1 
— these configurations have only graphs with tadpoles and bubbles on external lines, which 
vanish in dimensional regularization. 

In summary, for amplitudes with a single quark pair and idcntical-helicity gluon legs, 
there are two independent classes of primitive amplitudes to compute, 

A^-^{<P, 1/, 2+, . . . , j+, {j + 1)+, . . . , n+) and A'M 1/", 2+, • • • , j/, (j + 1)+, . . . , n+) . 
D. Spinor-helicity notation 

The primitive amplitudes are functions of the massive momentum A;^ of the particle, 
and the massless momenta ki of the n partons. All momenta entering the amplitudes are 
taken to be outgoing, so the kinematical constraints are, 

n 

k^ + J2ki = 0, (2.29) 
1=1 

kl = ml, k1 = 0. (2.30) 

The amplitudes are best described using spinor inner-products [72, 78] composed of the 
partonic momenta, 

(j = 0-|/+) = u_{k,)u^{ki) , \jl] = (j+r) = u^{k,)u_{ki) , (2.31) 
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where u±(k) is a massless Weyl spinor with momentum k and positive or negative chirahty. 
We use the convention [j I] = sign(A;°A;f) {I j)* standard in the QCD hterature, so that, 

{ij)\ji]^2ki-kj = Sij. (2.32) 

We denote the sums of cychcly-consecutive external momenta by 

K.^ ^ K + fc^Vi + • • • + ^i-i + k';. (2.33) 

where all indices are mod n for an n-parton amplitude. The invariant mass of this vector is 
Si...j = Kl„j. Special cases include the two- and three-particle invariant masses, which are 
denoted by 

Sij = Kj^ j = (/cj -|- A)j) — 2ki ■ kj, Sijk = (A)j -I- kj -\- kk) ■ (2.34) 

In color-ordered amplitudes only invariants with cyclicly-consecutive arguments need appear, 
e.g. and Sj^j_|_i^j+2. It is convenient to introduce the same combinations as eq. (2.33) 

but with the momentum added as well, 

with invariant-mass s<^,i...j = K^ -,„y Longer spinor strings will also appear, such as 

j 

{k-\^i...j\l-) = ^(A;m)[mi] , (2.36) 

m=i 

{k-\$i,...j,$i,...j,\l'') ^ 5] (^mi)[mim2](m2 ■ (2.37) 

For small values of n, such strings can be written out explicitly as 

{k-\{a + h)\l-) = {ka)[al] + {kb)[bl] , 
{k+\{a + b)\l+) = [ka]{al) + [kb]{bl) , 
{k-\{a + b){c + d)\l+) = {ka){a+\{c + d)\l+) + {kb){b+\{c + d)\l+) . (2.38) 

III. REVIEW OF RECURSION RELATIONS AND FACTORIZATION 
A. On-shell recursion relations 



Here we briefly review the on-shell recursion relations found and proved in refs. [36, 37]. 
For further details we refer to these papers. The on-shell recursion relations rely on general 
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properties of complex functions as well as factorization properties of scattering amplitudes. 
The proof [37] of the relations relies on a parameter-dependent "[/c, /)" shift of two of the 
external massless spinors, k and Z, in an n-point process, 

[k, I) : Afc ^ Afc - zXi , Xi ^ Xi + zXk , (3.1) 

where z is a complex number. The corresponding momenta (labeled by Pi instead of ki in 
this section) are shifted as well, 

pf pf{z)=p'; + '^{k-\r\r), (3.2) 

so that they remain massless, pl{z) — — pf{z), and overall momentum conservation is 
maintained. 

An on-shell amplitude containing the momenta pk and pi then becomes parameter- 
dependent as well, 

A{z) = A{pi, . . . ,pk{z),pk+i, . . . ,pi{z),pi+i, . . . ,pn). (3.3) 

When ^4 is a tree amplitude or finite one-loop amplitude, A{z) is a rational function of z. 
The physical amplitude is given by A{0). 

If A{z) — > as 2; ^ oo, as in suitable cases in tree-level QCD [36, 37, 38, 40], then the 
contour integral around the circle C at infinity vanishes, 

^£f^W^O; (3.4) 

that is, there is no "surface term". In appendix B we show that a tree amplitude A{z) 
containing an extra field also vanishes as 2; —> 00, for the same choices of shift that lead 
to a vanishing A{z) in pure QCD. 

Evaluating the integral (3.4) as a sum of residues, we can then solve for ^4(0) to obtain, 

A(0) = - V Res^. (3.5) 

^ — Z 

poles a 

If A{z) only has simple poles, then each residue is given by factorizing the shifted amplitude 
on the appropriate pole in momentum invariants [37], so that at tree level, 

A{0) = ^ A^(^ = Zrs)-^A-^\Z = Zrs) , (3.6) 
r,s,h '■•■■^ 
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where h = ±1 labels the helicity of the intermediate state. There is generically a double 
sum over momentum poles, labeled by leg indices r, s, with legs k and / always appearing 
on opposite sides of the pole. By definition, the leg I is contained in the set {r, . . . ,s} 
and the leg k is not. The squared momentum associated with the pole, K^...s, is evaluated 
in the unshifted kinematics; whereas the on-shell amplitudes Al and Ar are evaluated in 
kinematics that have been shifted by eq. (3.1) with z — Zrs, where 

To extend the approach to one loop [52], the sum (3.6) should also be taken over the two 
ways of assigning the loop to A^ and Aj^. This formula assumes that there are no additional 
poles present in the amplitude other than the standard poles for real momenta. At tree level 
it is possible to demonstrate the absence of additional poles, but at loop level it is not true. 

Because of the general structure of multiparticle factorization [79] , only standard single 
poles in z arise from multiparticle channels, even at one loop. However, double poles in z 
do arise at one loop due to coUinear factorization [52, 56]. The splitting amplitudes with 

helicity configuration (+ + +) and ( ) (in an all-outgoing helicity convention) can 

lead to double poles in z, because their dependence on the spinor products takes the form 

[ab] I {ab)^ for (+ + -1-), or its complex conjugate (ab) / [abf for ( ) [64]. As discussed 

in ref. [52], this behavior alters the form of the recursion relation in an essential manner. In 
general, underneath the double pole sits an object of the form, 

{ab) 

which we call an "unreal pole" because there is no pole present when real momenta are used; 
it only appears, as a single pole, when we continue to complex momenta. As we shall discuss 
in section VII, the finite (pqqgg ■ ■ ■ g amplitudes exhibit similar phenomena, except that in 
this case, just as in the pure-QCD case of finite qqgg ■ ■ -g amplitudes [53], one encounters 
unreal poles that do not sit underneath a double pole. 



B. One-loop factorization properties 

In order to build on-shell recursion relations, we need the factorization properties of one- 
loop amplitudes for complex momenta. It is useful to first review the factorization properties 
for real momenta, which we know from general arguments [72, 79, 80]. 
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As the real momenta of two color- adjacent external partons a and b become coUinear, the 
one-loop amplitude factorizes as, 

^(,) ^ ^ (SplitL°i(a'^% b'^; z) ••,(« + b)\ • • •) 

h=± 

+ SplitL^i(a'^«, z) A^U. . . , (a + b)\ ...))■ (3.9) 

The tree and loop sphtting amphtudes Spht*^°^ and Spht^^-* are given in ref . [64] . In general, 
there are three types of limits: a pair of gluons becoming coUinear; a gluon becoming coUinear 
with a quark (or anti-quark) ; and a quark and anti-quark becoming coUinear (if they happen 
to be adjacent). When all external gluons have positive helicity, as here, the limits simplify. 
For 0-amplitudes, no singularity results when a parton becomes coUinear with the state, 
because the (p is massive. 

When two gluons become coUinear, the finite (pqqgg ■ ■ ■ g amplitudes A^ and A^ behave 
respectively as, 

^^(0, . . . , a\ b\ . . .) SphtL") (a+, 6+; z) At,i<P, ••.,(« + 6)+, . . .) 

+ Split2^^^(a+, 6+; z) A^^U<I>, . . . , (a + 6)", . . .) , (3.10) 
Al{<f>, . . . , a+, 6+, . . .) Split (a+, 6+; z) AUi<P, ■■■,(« + &) + , ■ ■ •) 

+ Split«^^(a+, 6+; z) <°li(0, . . . , (a + 6)^ . . .) . (3.11) 

Here Split^^^^^ (Split^^)^^) is the one-loop splitting amplitude with a gluon (scalar) circulating 
in the loop. The remaining two terms, with opposite intermediate-gluon helicity, vanish 
because Split+'*(a"'", 6"*"; z) and A|j°2]^(0, ±, -|-, -|-, . . . , -|-) are zero. Because the two one-loop 
splitting amplitudes are equal [64], 

SplitJ^^^(a+, 6+; z) = SplitV'^"(a+, 6+; z), (3.12) 

we are motivated to take the difference between A^ and A^ to form A^~^. This combination 
has the simpler coUinear behavior, 

A^-(0, . . . , a+, 6+, . . .) ^ SplitL°V+, ^) An-Act>. . . . , (a + 6)+, . . .) . (3.13) 

When a gluon becomes coUinear with any fermion, both A^ and A* behave as (using 
fermion helicity-conservation and eq. (4.21)), 

A^'^(0, af, b+, . . .) Split|;5^(a±, 6+; z) A^1^(0, (a + b)f, . . .) . (3.14) 
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Hence the difference * lias tfie limiting behavior, 

af, &+, . . .) Splitj;5^(a±, b+; z) (a + h)f, . . .) . (3.15) 

When a quark and anti-quark become coUinear (for jf — 2), the amplitude factorizes onto 
the finite one-loop (f)gg . . . g amplitudes, 

A^'^(0, af, 6/, . . .) ^ SplitL°V7, (a + 6)+, . . .) 

+ Sphtf (a/, 6+; z) aK(0, (a + 6)", . . .) • (3.16) 

As mentioned in section II C, in contrast with the pure-QCD case [53], the gluon and scalar 
loop contributions to A„_i;i(0, ±, -|-, -|-, . . . , -|-) are not the same. For the difference A^~^ we 
find from eq. (3.16) the coUinear behavior, 

^^^(0, a7, 6+, . . .) ^ SplitL°V7, bl; z) AfX^. (a + 6)+, . . .) 

+ Sphtf (a/, 6/; z) A^^^cj,, (a + 6)", . . .) . (3.17) 

Next consider multiparticle factorization. In this case, the vanishing of the tree ampli- 
tudes An\lj ,2^ , . . . jjj}', . . . ,n^) and An\(f), 1^,2+, . . . . . . ,n~^) means that the only 

singularities are on gluon poles. There are two possible ways to attach the 4> field, 

A^'^{<P, 17, 2+, . . . , j/, . . . , m+, . . . , n+) ""-^^ (3.18) 
^^^,(0, 17, 2+,...,j;,..., m+, -K^_J-^Ail^^,{Kt.m, {m + 1)+, ■ ■ ■ , n+) 

-"■l...m 

+ A^Zii^^^ 2+, . . . , j/, . . . , m\ -K^ j'a^^:!1^,{<P, Kt„^, (m + 1)+, . . . , n+) , 

-"■l...m 

with m > 3 and m > jf. For A^~^, the first term in eq. (3.18) does not contribute, but the 
second does. 

In addition to the factorization properties for real momenta just described, the appearance 
of unreal poles of the form (3.8) sometimes has to be taken into account. While such behavior 
is not understood in general, in the present case we can use the similarity of our amplitudes 
to the corresponding finite pure-QCD amplitudes qqgg . . . g [53] . The contribution of unreal- 
pole terms to the recursion relations for those amplitudes could be stringently cross checked 
against various QCD and QED amplitudes. The inclusion of a massive (j) field is quite mild 
from the point of view of factorization, so it should not alter the form of the unreal-pole 
terms in the recursion relations. In section VII we will describe these terms in more detail. 
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IV. REVIEW OF KNOWN FINITE AMPLITUDES 



In this section, we collect previously-known results for tree and one-loop finite amplitudes 
in pure QCD, and for tree amplitudes containing a single (j) field. These results feed into 
the recursive formulae for the finite one-loop amplitudes containing a single (j) field, to be 
discussed in section VII. 

A. Pure-QCD amplitudes 

The n-gluon tree amplitudes with zero or one negative-helicity gluon all vanish, 

4'^)(1^2^3+,...,r^+) = 0, (4.1) 

where the omitted labels ( ". . ." ) always refer to positive-hehcity gluons. Also vanishing are 
the tree amplitudes for a pair of massless quarks and {n — 2) gluons all of the same helicity, 

4°)(17, 2+, . . . , (j - 1)+, j;, {j + 1)+, . . . , n+) = , (4.2) 

where the subscript / denotes a generic fermion, and the positive-helicity fermion (j^) can 
be located at an arbitrary position with respect to the negative-hehcity one (1^)- (Note that 
only the case j = 2 is required in the pure-QCD version of the tree- level formula (2.10).) 
The first nonvanishing tree amplitudes are the MHV amplitudes [72, 81], 

^<^*X,...,^,...,r^,...,r^)- . (4.3) 

and 

4°'(17,2^. .,m-,...,n-)^.^j|g^. (4.4) 

We also need the one-loop pure-gluon amplitudes with zero or one negative helicity. The 
all-positive-helicity case, for n > 4, is given by [59, 60] 

^|fVl+,2+,...,n+) =^W(l+,2+,...,n+) = - , , , , ^ — r, (4.5) 

' ' ' ' ' ' ' ' 3(12)(2 3)---((n-l)n)(nl) ' ^ ' 



where 



Hn — ^ ] Tr_ , (4-6) 

I<il<i2<«3<i4<i 
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and 

Tr. 



^ {1112) [12 {hi a) [i All] ■ (4.7) 



It is convenient to rewrite eq. (4.5) using an identity proved in ref. [53], 

1 



//„ = -^Tr_[^(2,r^)]^ (4.8) 

where we define the more general quantity T{l,p), 

p-i p 

i=l m=i+l 

involving sums over ordered pairs (i < m) of momenta between I and p. In eq. (4.8) for 
the n-point all- plus amplitude, all such pairs appear except those involving leg 1. Using the 
cyclic symmetry of this amplitude, any one of the n legs could be omitted, say leg Z, if the 
ordering in the corresponding sums for !F{1 + 1, Z — 1) is defined modulo n in eq. (4.9). 

For n = 3 we need to take one leg off shell. We define a vertex, with legs 1 and 2 the 
on-shell external legs [52], 



4^^(l+,2+,3+) = ^A^^^, (4.10) 

-^1,2 

where 

V3«(l+,2+,3+)^-^[12][2 3][31] . (4.11) 

The remaining one-loop finite helicity amplitudes in pure QCD are those for n gluons, 
with one of negative helicity; and those for a quark-antiquark pair plus {n — 2) gluons of 
positive hehcity. They are not actually needed for the recursion relations for the single-0 
amplitudes considered here. However, their structure, in the form given in ref. [53], is so 
similar to that of the single-0 amplitudes, that we reproduce them here. The one-minus 
n-gluon amplitude is given by, 

>l„.i(l-,2+,3+,...,n+) = — r, (4.12) 
' ' ' ' ' 3(12)(2 3)---(nl) ' ^ ' 

where 

^ _ :^ (io (i(/ + i))(i-|Am^(^+i)-Ji^) 
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n— 2 n— 1 



((/-i)0 



To = V V 



X 



(p 



(11 ^2...(/-l)^/...p |(P + 1)" 

x(l-|^,..,^(,+i)...„|l+>' 

(l1^2...(^-l)[-^(^>P)]'^(p+l)-n|l+) 



X 



The finite one-loop QCD amplitudes with quarks are [53], 



A^^(l7,2+,...,j/,...,n+) 



where 



and 



where 



2(12) (2 3) ■■■(nl) 



n— 1 



A^(l7,2+,...,j/,...,n+) 



S1 + S2 



Si 



So 



3(12) (2 3) •••(nl) 
ijl) (l(i + l))(l-|A^+i^a+i)-n|l+) 



n-l 

n—2 n—1 

{l-\tip^l)...ntl...p\il - (11 ^(p+l)...n^^..p 1^-^) 



((Z - 1) 



X 



{p (p + 1)) 



{1-\I/C2-(1-I)^l-P (11 ^2...(Z-l)^/...p I(P + 1) + ) 

X{^'\tl-p${p+l)-n\^^f {r\ $l...p$(p+l)...n |l^> 
(l-|^2...a-l)[^(/,p)]'^(p+l)...n|l + ) 



X 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



These expressions coincide for n = 4, 5 with previously-computed amplitudes [74, 82, 83]. 
Equation (4.12) agrees numerically with the result of ref. [60] for n < 18. Furthermore, 
in ref. [53] the amplitudes (4.15) and (4.17) were related to QED and mixed QED/QCD 
amplitudes by converting gluons into photons using appropriate permutation sums [74]. 
The resulting amplitudes can be compared to the one-loop QED and mixed photon-gluon 
amplitudes computed by Mahlon [60, 84], which provides a strong cross check of eqs. (4.15) 
and (4.17), and of the recursion relations used to generate them. 
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B. Amplitudes with a 4> field 



Amplitudes containing a single field have an "MHV structure" very close to that of 
pure QCD [23, 24]. First of all, those helicity configurations with almost all positive-helicity 
gluons, which vanish in tree-level QCD, eqs. (4.1) and (4.2), continue to vanish when a single 
(f) is added, 

4o)(0,l±,2+,3+,...,n+) = 0, (4.20) 
4°)(0,l7,2+,...,(j-l)+,j;,(j + l)+,...,n+) = 0. (4.21) 

Secondly, the amplitudes with a single field and two negative-helicity gluons, or a negative- 
helicity gluon plus a qq pair, and the rest positive-helicity gluons, are identical in form to 
the pure-QCD amplitudes, eqs. (4.3) and (4.4) [23, 24], 

>lf (0, 1+, 2+, . . . , mr, . . . , m,-, . . . , n+) = ^ 2)^73^ • > 1) ' ^^"^^^ 

and 

(Compared with the conventions in refs. [23, 24], the partial amplitudes here have an extra 
overall factor of i.) 

Equation (4.22) was asserted in ref. [23] based on comparison with exphcit results for 
small values of n [18, 22, 85], and consistency with collinear and multiparticle factorization. 
It can easily be proved recursively, however, by adapting the proof given in ref. [36] for the 
corresponding pure-QCD MHV amplitude. The pure-QCD recursion relation (3.6) is still 
valid for the amplitudes with a (p attached, because these amplitudes continue to vanish as 
2; — > oo, as we show in appendix B. 

For convenience, we take the negative-helicity gluons to be legs 1 and j. We use a [l,n) 
complex momentum shift; that is, we set A; = 1 and / = n in eq. (3.1), 

[1, n) : Ai ^ Ai - zXn , A„ ^ A„ + zXi . (4.24) 

Then the on-shell recursion relation [36, 37] has only one nonvanishing term, depicted in 
fig. 3, 

Af(0,l-,2+,...,r,---,ri+) 

= 4°1,(0, i-, 2+, . . . , . . . , (n - 2)+, i:+_i,n) Af\{n - 1)+, n+, -i^-_,,J . 

^n— l,n 

(4.25) 
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FIG. 3: Diagram corresponding to the recursion relation (4.25) for An\cl), 1 , 2+, . . . , j , . . . , W^). 

The form of this relation is precisely the same as for the MHV ra-gluon amplitudes in pure 
QCD: The only place the shift has any effect is in Kri-i,n, and it has no dependence on 
the (implicit) momentum of 0. Therefore, the solution of the recursion relation is exactly 
the same as in pure QCD, for n > 3. We demonstrate the solution (4.22) recursively by 
substituting it into the right-hand side of eq. (4.25), and recovering the left-hand side, 

"4 3 

(Ij) i [{n-l)n] 



(1 2) (2 3) ■ ■ • ((n - 2) K^_,,r.) (^n-l.n 1) ^n-l,n [u (-ir„_i,„)] [(-i^n-l,™) (u - 1)] 

' (1 2) ■ ■ ■ ((n - 3) (n - 2)) {{n - 1) n) (l"! |(n - 1)-) {{n - 2)-| l^') 



(4.26) 



(12) (23)---(nl) 

There is one important exception at the lowest order, n = 3. Here the recursion for the </> 
case has one extra term, replacing the term in eq. (4.25), which vanishes for n = 3, 

4°)(</),l-,2-,3+) = Af\cP,i-,k^^,)—AP{2-,3\~kt,3)- (4.27) 

This term is nonvanishing, and the three-parton kinematics are nonsingular, because of the 
extra momentum injected by 0. The proof of the fermionic MHV formula, eq. (4.23), is 
completely analogous. 

Another infinite series of tree- level (p amplitudes is the "maximally googly" set [24], 

4°'(0.1-.2-.3-.....«-)=,^gJ^. (4.28) 

A^(,U\2^X,...,n*)-^j^^j^g^y (4.29) 
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C. Coupling of A to the axial current divergence, and the (pg^g'^ amplitude 

As mentioned in section II A, at order in the effective theory, the pseudoscalar A 
begins to couple to the divergence of the hght-quark axial current [69]. The full effective 
Lagrangian is given by 



where the sum is over all light quark flavors rii, and C and C denote the appropriate 
coefficients, which begin at order for C and order for C. It was found in ref. [69] 
by explicit calculation that the second term in eq. (4.30) does not contribute at NLO, in 
accordance with a non-Abelian version [68] of the Adler-Bardeen theorem [67]. This can 
also be seen via recursion relations as follows: The amphtude Af\Aj^, IJ, 2^) vanishes by 
helicity conservation for massless fermions /. We have denoted the contribution from the 
second term in eq. (4.30) by the subscript J5. By A^^ we mean a tree amplitude, which 
contributes at NLO due to the factor of in C The effective Agg- vertex is proportional to 
the momentum that flows through it due to the derivative in eq. (4.30). When attaching a 
gluon to this vertex, one can choose its polarization vector so that each of the two diagrams 
contributing to Af\Aj^, IJ, 2^, 3^) vanishes. (The amplitude is of course gauge invariant; 
in other gauges the two diagrams cancel against each other.) Thus Af\Aj^, IJ, 2^ , 3^) = 0. 
By recursion, the amplitudes with an arbitrary number of gluons also vanish — there are 
no lower-point amplitudes to contribute to the right-hand side of a recursion relation. Thus 
we do not need to consider the contribution of the light-quark axial current at NLO, and 
the amplitudes involving a pseudoscalar A can be constructed via eq. (2.7). 

In appendix A we compute the simplest, finite one-loop amplitude with a single (p field, 
^2;i(0, 1"*", 2''"). The result can be written in terms of the "googly" two-gluon tree amplitude. 



This amplitude enters recursively into all the subsequently constructed amplitudes. 
V. THE SOFT HIGGS LIMIT 




(4.30) 



1=1 



^2;i(0,l+,2+) = -2^^'^^(0t,i+,2+). 



(4.31) 



It is useful to consider the kinematic limit of the single-0 amplitudes as /c^ 0. Of course 
this soft-Higgs limit can only be reached if uih — 0. The limiting behavior of single-0 tree 
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amplitudes was studied in ref. [23]. Here we review that analysis and extend it to the loop 
level. First consider amplitudes with a single real scalar field H , as the Higgs momentum 
goes to zero, so that H has almost no spatial variation. Then the HiiG^,^G^'' interaction 
behaves just like a constant multiple of the pure-glue Lagrangian trGfj^^G^". So one might 
expect the Higgs-plus-n-gluon amplitudes to become proportional to the pure-gauge-theory 
amplitudes, according to the relation, 

d 

An{H, {ki, Xi, Qi}) — > (const.) X —AnHh, Xi, Qi}) askH^O, (5.1) 

dg 

for any helicity configuration. Taking into account the number of factors of the gauge 
coupling in the n-parton tree amplitudes, namely (n — 2), eq. (5.1) leads to the tree- level 
soft Higgs hmit, 

A^^\H, 1,2,3, .. .,n) — > (n - 2)A(f)(l,2,3,...,n) askn^O, (5.2) 

which holds color structure by color structure. Similarly, at I loops, our naive expectation 
for the soft limit is, 

A^')(i/,l,2,3,...,n) — ^ (n + 2i-2)^(,')(l,2,3,...,n) as kn ^ 0- (5.3) 

In fact, wc will find that eq. (5.3) is not precisely correct beyond tree level. A generalization 
of this equation to the case of 4> and (p^ amplitudes will hold at loop level, but again it will 
hold precisely only for the finite one- loop amplitudes. Equation (5.3) generically breaks down 
at loop level because of an exchange-of-limits problem, between dimensional regularization 
of infrared divergences, and the momentum injected into the amplitude by the Higgs boson, 
which acts as a kind of effective infrared regulator as well. 

Before examining the infrared issue, let us first propose an extension of eq. (5.3) to the 
case of amplitudes with a single 4>, or 4>^. At tree level, by studying the MHV representation 
of the (f) and 0^ amphtudes, it can be shown [23] that the factor of {n — 2) in eq. (5.2) is 
partitioned into (n_ — 1) and (n+ — 1) in the two cases, 

A^°^(0, 1,2,3, ...,n) — > (n_ - 1)4°) (1,2, 3,..., n) ask^^O, (5.4) 

and (by parity) 

4°)(</.t, 1, 2, 3, . . . , n) ^ (n+ - l)Af (1, 2, 3, . . . , n) as A;^ ^ 0, (5.5) 
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FIG. 4: Two cuts of a one-loop (p amplitude which become the same cut in the soft limit k(j) — > 0: 
(a) the cut in the Si...j channel, and (b) the cut in the S(f,^i...j channel. 

where n± are the number of external partons with helicity ±1 for gluons (±| for fermions). 
The natural generalization of eqs. (5.4) and (5.5) to / loops would be, 

All\(j), 1,2,3, ... ,n) — > (n_ + / - l)4')(l,2,3,...,n) ask^^O, (5.6) 

and 

A^J^\(f)\ 1,2,3, ... ,n) — > (n+ + /-l)A«(l,2,3,...,n) ask^^O. (5.7) 

But again, infrared complications can affect these relations. 

Consider a particular class of terms in a one-loop Higgs or (p amplitude, those containing 
a cut in either the Si...j channel or the s,j)^i...j channel, as illustrated in fig. 4. The soft behavior 
of these terms can be analyzed using unitarity, plus the tree-level soft behavior. The two 
cuts in fig. 4 become identical, up to constant factors, as k^ 0. Suppose there are m 
external partons on the right side of the cuts in fig. 4, plus the two crossing the cut. Then 
there are (n — m) external partons on the left side of the cut, plus two more crossing the 
cut. (The fields H and do not count as partons.) For the case of a Higgs field H, using 
eq. (5.2) in the soft limit kff 0, the cut in fig. 4(a) reduces to the pure-QCD cut in the 
Si...j channel, multiplied by a factor of {n — m). The cut in fig. 4(b) reduces to the same 
pure-QCD cut, because s^^i...j —>■ Si...j, but it is multiplied by a factor of m. Summing the 
two contributions, we find that the terms in the loop amplitude containing ln(— Sj...^) and 
ln(— S(/,,j...j) go smoothly into the terms containing ln(— Sj...j) in the pure-QCD amplitude, 
multiplied by a factor of {n — m) + m = n. This argument would seem to prove eq. (5.3) for 
/ = 1. 
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A similar argument applies to the and (f)^ amplitudes. If there are a total of n± x, 
partons with helicity ±1 on the left-hand side of the cut, and n±^R on the right-hand side, 
then n±^L + n±^ii — n± + 2, because 2 partons of each helicity are assigned to the cut. Instead 
of the multiplicative factor oi {n — m) + m — n found for the case of if, for we get from 
eq. (5.4) a factor of {u-^l — 1) + {i^-,R ~ ^) — i^-i which is the factor in eq. (5.6) for I — 1. 

However, both these arguments break down for the case m = 2, for which only two 
partons, i and i -|- 1, appear on the right-hand side of the cut. In this case, the cut in 
fig. 4(a) is generically infrared divergent, whereas the cut in fig. 4(b) is infrared finite, for 
nonzero momentum k^. The generic form of one- loop divergences for the unrenormalized 
amplitudes for an if or field plus n gluons, for example, is [86] 



^„;i(ii,l,2,...,n) = A^^\H,l,2,...,n) + finite, 

^„;i(0,l,2,...,n) = - 4°n0,l,2,...,n) + finite. 

The infrared divergences have precisely the same form as for pure QCD, 



(5.8) 
(5.9) 



>ln;l(l,2,...,n) 



A(f)(l,2,...,n) + finite, (5.10) 



where /3o = (H-^c — 2n/ — ns)/{SNc). (After ultraviolet renormalization the remaining 
infrared Pq terms match between eqs. (5.8)-(5.10).) Hence the divergent parts of the (renor- 
malized) one-loop H and amplitudes have the soft behavior characteristic of the tree-level 
amplitudes. 



1/ 1 ' ^^-^^^ 

l/e poles 



^„;i(ii,l,2,...,n) (n-2)A,.i(l,2,...,n) 

l/e poles 

A„.i((/), 1, 2, . . . , n) (n_ - 1) A„.i(l, 2, . . . , n) . (5.12) 

l/e poles l/e poles 

This behavior docs not agree with the naive expectation for the one-loop behavior, eqs. (5.3) 
and (5.6) for / = 1. However, it is consistent with the behavior of the divergent cuts for 
j = i -I- 1 (m = 2) in fig. 4. For example the factor of {n — 2) in eq. (5.11), for the cut of 
the term proportional to ln(— Sj^j+i)/e x An \ comes just from fig. 4(a), because the cut of 
fig. 4(b) has no pole in e. It is the exchange of hmits, e — > and kn — > 0, in fig. 4(b), which 
is at the root of the non-uniform soft behavior. 

This argument does not cover the cases m = 1 (ln(— j) terms) or m = (ln(— A;|) terms). 
These cases have only the cuts shown in fig. 4(b), but the cuts become quite singular in the 
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/c<^ — limit, so we may expect non- uniform soft behavior. The cut-based argument also 
does not tell us what to expect for the rational terms in generic H oi (j) amplitudes. Indeed, 
for the one-loop HqqQQ amplitudes presented in ref. [21], the soft Higgs limits of the generic 
cut terms (all of which have m — 0, 1, or 2) and rational parts at order e° do not appear to 
be uniform. 

On the other hand, for the finite amplitudes we focus on in this paper, there should 
be no infrared issues, and hence we should find the simple, uniform behavior, valid for each 
color structure, 

A^j^\(j), 1, 2, 3, . . . , n) ^'^ °> n_ ^^^■'(1, 2, 3, . . . , n) (finite helicity configurations). 

(5.13) 

We shall verify eq. (5.13) explicitly below for the finite helicity amplitudes. 

Assuming that eq. (5.13) holds, then the generic one-loop amplitude with (p present 
must be at least as complicated as the corresponding pure-QCD amplitude with omitted, 
because the former amplitude must reduce to the latter amplitude, multiplied by a factor of 
n_, as 0. The one possible exception is the all-plus amplitude, An\(t), 1^, 2+, . . . , n+), 
because the multiplicative factor vanishes in this case. Indeed, as we shall see in the next 
section, these amplitudes are given by a much simpler formula than the corresponding QCD 
result (4.5). 



VI. RECURSION RELATION FOR ALL-PLUS CASE 

Using the results of section 111, and the explicit calculation of the two-gluon one-loop 
amplitude (4.31), we can recursively construct the following form of the one-loop all-plus 
amplitude, 

4^1(0,1+, 2+,..., n+) = AW(0,l+,2+,...,n+)=O, (6.1) 
A[fl(0,l+,2+,...,n+) = -2z J^^^^ _^^^^ = -24°) (0t, 1+ 2+, . . . , n+) . (6.2) 

The proof of eq. (6.2) is analogous to the above proof of the tree-level MHV amplitude via 
an on-shell recursion relation [36]. 

Here we use again the [l,n) shift (4.24). As in section IV B there is only one term 
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FIG. 5: Diagram corresponding to the recursion relation (6.3) for An {<p,l~^ ,2~^ , . . . , n~^). 
contributing to the recursion, as illustrated in fig. 5, 

= A„_i;i(0, i+, 2+, . . . , (n - 2)+, K^,,J APiin - 1)+, n+, -K-_,^J ■ (6-3) 

It is straightforward to show that eq. (6.2) satisfies eq. (6.3), using the same algebra as 
m eq. (4.26). The amplitude aL^'(0, 1+, 2+, . . . ,n~^) in eq. (6.2) vanishes in the soft limit, 
just as predicted by eq. (5.13) for n_ = 0. 

Note that the lack of ny or Ug dependence in A2-i{(f),l^ ,2^) translates recursively into 
the vanishing of aK^ and An'^ in eq. (6.1), for all values of n. Now consider the dependence 
of the one- loop one- minus (p amplitude on nf or n^. This amplitude factorizes onto two 
types of one- loop amplitudes, either the all-plus amplitude, which obeys eq. (6.1), or the 
all-plus pure-QCD amplitude. But the pure-QCD amplitude obeys An^ = —An\ thanks 
to a supersymmetric Ward identity [76]. Using this relation and eq. (6.1) in a recursive 
argument, we must have An'^ = —An'^ also for the one-minus amplitude, as mentioned in 
section II, and as assumed in the color decomposition (2.18). 



VII. QUARK RECURSION RELATIONS 

As we discussed in section II, for the finite helicity amplitudes — one quark pair, with 
all gluons having positive helicity — there are two independent primitive amplitudes that 
need to be computed, 

A^-%<j),j;) ^ A^^(0,l7,2+,...,j/,...,n+), 2<j<n, (7.1) 
A:(0,j/) = A:,{<f),l-,2+,...,j;,...,n+), 2<j<n-2. (7.2) 
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We abbreviate the notation as in ref. [53], retaining only the label of the positive- helicity 
fermion (and 0). A computation of these primitive amplitudes determines all the finite 
one-loop quark amplitudes for a single (j) boson coupled to QCD. 



A. Recursion relation for L — s contribution 

As mentioned in the introduction and in section III, finite quark- gluon amplitudes in both 
pure QCD and (/)-amplitudes contain "unreal" poles. The contributions of the unreal poles 
can be described using a "soft factor" iS'-'^-* reminiscent of the propagation of an internal soft 
gluon [53]. 

For the recursive construction we need the "L — s loop vertex" introduced in ref. [53], 

F3^-^(l+,2+,3+)^-^[12][2 3][31] . (7.3) 

The soft factor S^^\a,s,h) traditionally describes the insertion of a soft gluon s between 
two hard partons a and 6 in a color-ordered amplitude. It depends only on the helicity of 
the soft gluon and is given by [72], 

S^'\a^s\h) = (7.4) 

(a s) {s 0) 

S^'Ha,s-,b) = -T^y (7.5) 

[a s\ [s b\ 

Choosing the shifted legs in eq. (3.1) to be [k, I) = [1, n), eq. (4.24), the recursion relation 
for the pure-QCD amplitude A^^^{j^) = A^~^{1^, 2+, . . . , j^, . . . , n'^) was found to be. 



= A^Ztiii, 2+, . . . , j/, . . . , (n - 2)+, X+_,_J — — 4'^(-X-_i,„, {n - 1)+, n+) 



Sn—l,n 

|(0) 



+ ^i"!,(i7, 2+, . . . , j/, . . . , (n - 2)+, K-_,^J Vt\-KU,n. (n - 1)+, n+) 

^n— l,n 

X 5(°)(j, i) S^'\n, n-1), (7.6) 

where the hatted momenta in eq. (7.6) are defined by the shift (4.24), with 

^l-i,n _ ((n-l)n) 
(l-|^„_i,Jn-) ((n-l)l)' 

in each term. The second term in cq. (7.6) represents the contribution of an unreal pole. 
The solution to eq. (7.6) is given by eq. (4.15). 
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FIG. 6: Diagrams corresponding to the terms in the recursion relation (7.8) for A^~^{(j),jj~), for 
j < n — 2. The second term contains an unreal pole. In the last term, the sum over I runs from 
j + 1 to n, inclusive. 

For the analogous smgle-0 amplitude 74^~'*(0, j^), using the same [1, n) shift, the n-point 
on-shell recursion relation is, for j < n — 2, 



A^Z^ (0, 17, 2+, . . . , j;, . . . , (n - 2)+, k^_,^J Af (n - 1)+, n+) 

+ 17, 2+, . . . , j/, . . . , (n - 2)+, K-_,J Vt\-Kt-i,n, {n - 1)+, n+) 



n-l 
n 

+ 5^ A^ (17, 2+, j/, ...,(/ - 1)+, i^7,.. J ^ 41+2(0, -^,\..n> n^) , 

(7.8) 

where in the s<^,;...„ channel we set. 

The recursion relation is shown in fig. 6. Again the second term represents the unreal-pole 
contribution, which we take to have the same form when is present as it does in pure 
QCD. The additional third term in eq. (7.8), containing a sum over /, arises because the 
amplitudes Am {(f), +,+,...,+) = Am~'^\(f), +,+,...,+) are nonvanishing when is present. 
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The corresponding pure-QCD amplitude combination Am ^\+, +!•••,+) vanishes according 
to eq. (4.5), so there is no such term in the pure-QCD recursion relation (7.6). 
The solution to eq. (7.8) is, 

where 

n-l 

Qt = (lJ)E(l1^'^^-|l^> ' (7-11) 
Qt = _4 f 

X ^2...a-l)^^...n |l + >' {r\^l-n^l-n |l + > , (7-12) 

Qt^-ATT ^ii^ 

l^^,tltl (11 ^p...n(^l...a-l) + K-n) \il - 1)^) (11 ^p...n(^l...(;-l) + ^p-n) 

^ ((P- 1)P) 

(1-|^2...(Z-1)(^1...(/-1) +^p...„) \{P-Ir) {l~\t2..il-10l-(l-l)+t,-n) 

Making use of momentum conservation, eq. (2.29), antisymmetry of the spinor products, 
and re- indexing the sum over p in Q2, we can rewrite these quantities as 

n-l 

Qt = (1 i)E (11^^,2...^' 1 1-")' (7-14) 

1=2 

,4^1 (i"l ^i-Mi^ - 1)^) n 

X |l+>' ^,.J^ , (7.15) 



^ {P (P+1)) 

^2...a-l)^0,Z...p (11 ^2...(Z-l)^0,Z...p I(P + 1) + ) 

X (1-| |l+>' (1-| |j+> ^ . (7.16) 

The latter forms make the soft-0 limit, A;^ — 0, somewhat more manifest. In this limit, we 
see that Qq reduces to Qo as given in eq. (4.16), and that Qf and Q2 both vanish. (Note 

35 



that there is one more factor of in the numerator than in the denominator of Qi ) Thus 
we have, 

^n-^(0, J?) A^-'ij/) ask^^O, (7.17) 

in agreement with eq. (5.13) for n_ = 1. 

The recursion relation (7.8), and thus its solution (7.10), is strictly valid only for j < n—2. 
However, the cases j = n — 1 and j = n can be extracted from amplitudes with a larger 
number of legs n and the same value of j, via appropriate coUinear limits. In the limit 
that gluons n — 1 and n become coUinear in the amplitude A^~^{(f), {n — 2)^), we obtain the 
amplitude A^li(0, (n - 2)^) via eq. (3.13). Similarly, taking the fermion n—1 and gluon 
n to be coUinear in the amplitude A^~^{(j), {n — 1)^) gives the amplitude A^r^(0, {n — 1)^) 
according to eq. (3.15). Taking these limits, we find that eq. (7.10) continues to be valid for 
j > n — 2. That is, Q2 vanishes for j = n — 1, and Qf and Q2 vanish for j = n, as can be 
seen from the limits of the sums. 



B. Verification of solution for ^{(f),j^) 

In this subsection we demonstrate that the set of amplitudes A^~^{(f),j^) satisfies the 
recursion relation (7.8). The analysis parallels section VD of ref. [53], which verified the 
solution A^j(j^) = A^^{lj^ , 2"*", . . . . . . , n+) in cq. (4.17) for the scalar-loop contributions 

to the pure-QCD finite quark amplitudes. However, the present analysis is simpler, because 
the all-plus gluon amplitudes (6.2) with a single (p present are simpler than those in pure 
QCD. The first term on the right-hand side of the recursion relation (7.8), shown in fig. 6, 
includes the three-point tree amplitude Af\—K~_i^^, (n— 1)+, fi+). Let Qf, i = 0, 1, 2, stand 
for the shifted versions of Qf for the appropriate {n — l)-point one-loop quark amplitude, 
A^Zl{(p, if, 2+, . . . . . . , (n — 2)+, „). The first term in eq. (7.8) can be simplified to 

_1 Qt + Qt + Qt ^_ 

2 1) (1 2) (2 3) • • ■ ((n - 3) {n - 2)) {{n - 2) kn-i,n) ^n-i,n 

\(n — 1) n]"^ 



[n Kn-l,n] [Kn-l,n [n - 1)] 

1 Qt + Qt + Qt {{n - 2) (n - 1)) {{n - 1) n) {n 1) [jn - 1) nf 

2 (1 2) (2 3) ■ ■ ■ {n 1) {{n - 1) n) [{n - 1) n] (l"! |(n - 1)") {{n - 2)-\$^_,, 



\n 



Qt + Qi + Qt X8) 



2(12) (2 3)---(nl) 
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Thus the correct spinor denominator factor is reproduced. 

Next wc need to determine how the quantities Qf arc affected by the shift (4.24). First 
consider Qq, as given in eq. (7.14). The {n — l)-point expression Qq is a single sum over 
I containing {n — 3) terms, which is one fewer than the number of terms in the n-point 
expression on the left-hand side of the recursion relation (7.8). All terms except the last 
in Qo behave simply under the shift (4.24) of Ai and A„. They depend on Ai through 
and |1+), but are independent of Ai. The dependence on A„ is solely via the factor 

(l-|(^(;+l)...(„_2) + in-l,n) ■ ■ ■ = (l~l^a+l)...n ■ ■ ■ , (7-19) 

because the shift in $^n-i,n is proportional to Ai. So each such term directly yields the 
corresponding term in the n-point sum Qq. 

The missing last term in Qg (with / = n—1) is provided by the unreal-pole term containing 
V^^^{—K^_i^j^, {n — l)+,n+) in the recursion relation (7.8). This term can be simplified to, 

, 2 



i (1 Kn-l,n) {j Kn-l,n) [Kn-l,n (n " 1)] [{n - 1) u] [u Kn-l,n\ 



(1 2) (2 3) • • • {{n - 3) (n - 2)) {{n - 2) Sn-i,n 
(j 1) [n (n - 1)] 

{jkn-l,n) {kn-l,nl) (n - 1)] 



(liXil 


kn-l,n 


\n ) [{n — 1) n] 


2(12) (2 3)---((n-3) (n - 


-2)) {{n-l)n){{n-2)-\ 


kn-l,n 


\n-) 



i {lj){ln) [n (n-l)]((n-l) 1) 
2 (12)(2 3)---(nl) 

i(lj)(l-|^^,2...(n-l)4-l|l+) 



(7.20) 



2 (12)(2 3)---(nl) 

which clearly agrees with the term with I — n — 1 in the expression (7.14) for Qq. 

Note that the same algebra, with set to zero, demonstrates that the purc-QCD am- 
plitude A^-'{jp in eq. (4.15) solves the corresponding recursion relation (7.6). Com- 
pared to the single-^ recursion relation (7.8), the pure-QCD relation (7.6) lacks the 
set of terms containing the amplitudes for a particle plus all positive-helicity gluons, 
A^n-i+2(^i ~-^^i-ni ^'^f ■ ■■> ^^)- Thus these latter terms should serve as "sources" for the Qf 
and Q2 parts of the A^'~^{(f), jj') solution. 

Let us now turn to Q2, as given in eq. (7.16). Wc split the terms in Q2 into those with 
p < n — 1 and those with p = n — 1. The terms with p < n—1 come from the (n — l)-point 
Q2 contribution in eq. (7.18). To show this, we observe that, just as for the Qq terms, Ai 
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never appears in Q2. Also, A;„ only appears in eq. (7.16) for via . . .. Thus 

we may use eq. (7.19) once again (with / replaced by p), in order to see that every term with 
p < n — 1 in is generated directly from the corresponding term in Q2. Similarly, every 
term in with I < n comes from the corresponding term in Q^. 

The Q2 terms with p — n — 1, and the Qf term with I — n, originate instead 
from the terms containing the amplitudes for a particle plus all positive-helicity gluons, 
A|^1^^2(0) ~^4,i---ni 1 ■ ■ ■ i"^^)) in 6q. (7.8). The term in Q2 comes from the /*'^ term in 
the recursion relation, for / < n. For / = n, it gives the / = n term in . 

To write the t"^ term in the proper form, we use the following identities, for I < n, 

(nl) S^^i...n = (l~| ^2...(/-l)^<^,i...(n-l) , (7.21) 



(1 I $ct>,l-n 



= -•S(/.,Z-(n-l) , (7.22) 



{{n-l)n) = ((n-l)n)- ^ _. ((n - 1) 1) 

^ (l-|^2...(Z-l)^0,Z...(n-l)|(n-l) + ) 
(l"l^0,i-nl^") 

Then the V'^ term in the recursion relation (7.8) becomes, 

2 

{lK^^l...n) {jK^^i...n) 1 



(7.23) 



2i. 



(1 2) (2 3) • • • ((/ - 2) (/ - 1)) ((/ - 1) K^,i...r,) s<p,i...n 



X 



{K^,i...r^ /) (/ (/ + 1)) ■ ■ ■ ((n - 1) h) {n K^j...^ 
2i {{I - 1) {{n - 1) n) (n 1) (1-| \n-f {r\t4>,i-n 



n 



(1 2) (2 3) • • • (n 1) (n+l |t^,,.., |(/ - 1)+) |/+) 

(nl) (l-|^^,,..Jn-) (/c2)2 



X 



(1 I $2-{l-\)t<j>,l...{n-l) W) (1 I t2-{l-l)t4>,l-{n-l) I (^ - 1) + ) S</.,Z-(n-l) 
2i ((/-I)/) 

(1 2) (2 3) • • • (n 1) (1- 1 ht4>,i-in-i) I (/ - 1)+) (1- 1 Ltc^,l..in-l) 

((n — 1) n) 



X 



(1 I t2-{l-X)t4,.l...{n-\) l(n - 1) + ) (1 I t2-{l-\)t^,l...(n-\) W) 

X |l+>' 4^<^,^..(n-l) |j+> -^^^ • (7-24) 

-50,i-(n-l) 

The final form is just the term with p = n — 1 in formula (7.16) for Q\. 
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The case I — n is a little different. The relevant identities (7.21)-(7.23) are slightly 
modified, to 



(nl) s^^„ = (l I ^2-(n-i)l^<^ , (7.25) 

= -kl . (7.26) 

We simplify the Z*^ term in the recursion relation (7.8), for Z = n, as follows, 

_2 . {lK4>,nf{jK^,n) 1_ {kl? 

' (1 2) (2 3) • • • ((n - 2) (n - 1)) ((n - 1) S0,n {fi K^,r^f 

^ 2i {{n - 1) n) {n 1) (I'l ji:^ \n-)^ Q- 1 1^-) 
(12)(2 3)---(nl) (n+||^J(n-l)+) 



X 



(l-|^2...(n-l)^<^|ri+) {kir 

2/ 



(7.27) 



(1 2) (2 3) • • . {n 1) (1-| |(n - 1)+) (l"] 1^^) ' 

The final form is identical to the term with / = n in formula (7.15) for Q'^. We have now 
accounted for all the terms in + Q2 that were not present in Qi + Q2 > thus completing 
the proof that eq. (7.10) obeys the recursion relation (7.8). 

C. Recursion relation and solution for 

The recursion relation for the pure-QCD primitive amplitude ^^(jj") reads [53], 
^:_,(i7, 2+, . . . , j;, . . . , (n - 2)+, J Af (-it-_i,„, (n - 1)+, n+) 

*n— l.n 



+ Af_,{ij, 2+, . . . , j;, . . . , (n - 2)+, k-_,^^) V,^'^ {-KLi,n, (n - l)\ n+) 

^n-l,n 

X (1 + sn-i,n S^'^ (n - 2, j) (n, n - 1)) 



n-2 



+ E 4°^(i/^2^...,j;,...,(i-l)^^,-J^4V,(-i^+„,i+,...,n^ 

l=J+l 



(7.28) 
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FIG. 7: Diagrams corresponding to the terms in the recursion relation (7.30) for A^{(l),j^). The 
second diagram contains a double pole, and an unreal pole underneath it. In the last diagram, the 
sum over / runs from j + 1 to n — 2, inclusive. 

The hatted legs undergo the shift in eq. (4.24), where in the channel we set z to the 
value 

(7.29) 



(1 \^i-n \n ) ' 

The recursion relation for the corresponding single-0 amplitudes, A^(0,jj"), depicted in 
fig. 7, is very similar. 



^ ,(0) 



^:_,(0, 17, 2+, . . . , j;, . . . , (n - 2)+, ir+_,,J {n - 1)+, n+) 



+ A':U<P, 17, 2+, . . . , j/, . . . , (n - 2)+, K-.,,J V^'>{-Kti,n, (n - 1)+, n+) 

(1 + 5 W(n - 2, Kti,n,j) -R-.,^^, n-l] 



X 

+ E V' 2+, j;, ...,(/ - 1)+, i^,:. J Ai\,i-kt.^, n+) . 

(7.30) 

Notice that in this case only makes its appearance in the "source terms" — the second and 
third terms in fig. 7 — through tree amplitudes. Recall also that the relevant 0-containing 
MHV tree amplitudes (4.23) are identical in form to the ones without 0, eq. (4.4). Thus we 
should expect the solution to the recursion relation to have essentially the same form as in 
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the pure-QCD case. 

Indeed, the solution to eq. (7.30) is 



where 



St 



X 



^^ (l1^^..n[^^(^,p)]^^(p+l)■■■n|l+) ^^33^ 

The S'f term manifestly has the same form as the pure-QCD result in eq. (4.18), 
apart from the implicit (j) momentum. For the 5*2 term, we use momentum conservation to 
substitute 

^2...(/-l) = -ti...n , (7.34) 

in the expression for 5*2 in eq. (4.19). Then the 5*1 term is identical in form to 5*2 as well. 

Because of the simple relation between Sf and Si, it is trivial to show that the so- 
lution (7.31) for A^{(f),j^), like the solution (7.10) for A^~^{(f),jj'), obeys the expected 
soft-Higgs hmit (5.13) with n_ = 1, 

K{<f>,j;) At^U;) as ^ 0. (7.35) 

The simple relation between 5*^ and Si also makes it very easy to demonstrate that the 
set of amplitudes A^{(f),j^) given in eq. (7.31) satisfy the recursion relation (7.30). We just 
use the analysis in section VD of ref. [53], which proved the analogous result for QCD, 
namely that Af^{j^) as given in eq. (4.17) satisfies the recursion relation (7.28). The only 
modification that needs to be made to the arguments and identities is that ^2 - (/-i) should 
be replaced by — everywhere, according to eq. (7.34). After doing that, all expressions 
appearing in the two recursion relations and solutions, for A^ with and without, have 
exactly the same form. The momentum of (j) only appears implicitly in An{(j),jf~), through 
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momentum conservation. In addition, after replacing ijC^ 



2---(«-l)) 



none of the identities used in 



ref. [53] require momentum conservation (which could have introduced the momentum). 
In the analysis of the contribution of the V'^ term in the last diagram of fig. 7, for example, 
only legs I to n are involved in an essential way. 

D. Factorization properties of solution, and a (?!)-plus-n-gluon amplitude 

In appendix A of ref. [53] it was shown that the pure-QCD formula (4.17) has all the 
correct multiparticle poles, factorizing properly onto products of the quark- containing MHV 
tree amplitudes (4.4) and the one-loop all-plus pure-gluon amplitudes (4.5). The corre- 
sponding demonstration for the (j) amplitude A^{(l),j^) is completely analogous, so we do 
not present it here. The coUinear singularities of eqs. (7.10) and (7.31) work in exactly the 
same way; there are no coUinear singularities between the particle and another parton, 
only between pairs of partons, and the corresponding (n — l)-point amplitudes again have 
exactly the same form as in the pure-QCD case. By taking the anti-quark and quark to be 
coUinear in the amplitude in which they are adjacent (j = 2), we can extract the scalar and 
[g — s] loop contributions to the one-loop amplitudes for (pg'^g'^ . . .g^. 

A term containing 1/ [1 2], corresponding to the splitting amplitude Split?'* (1 J, 2^ ; z), in- 
dicates a factorization onto the all-plus amplitude, (f)g'^g^ ■ ■ ■ g'^. We observe that eq. (7.31) 
for ^4^(0, 2^) does not contain a factor of 1/ [12], thus confirming the observation in sec- 
tion VI that the scalar 0-amplitude with gluons of positive helicity vanishes, eq. (6.1). 
Similarly, in this coUinear limit, A^~*(0,2y) as given in eq. (7.10) factorizes onto the glu- 
onic 0-all-plus amplitudes, eq. (6.2). The relevant contribution comes from the term with 
I — 3 and p — n oi Q2, as given in eq. (7.13). 

Analogously, we can extract the scalar and [g — s] contributions to the amplitude for a 
plus n gluons, one of which has negative helicity, from eqs. (7.31) and (7.10), respectively. 
The result for the scalar contribution is. 



AW(0,1-,2+,3 



+ 




i Tf + 
3(12) (2 3)---(nl) 



(7.36) 



where 



n-l 



(iO(i(/ + i))(i-|^,,,+,^(,^i)...ji+) 

{I {I + 1)) 



(7.37) 



1=2 
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n— 2 n— 1 



({1-1)1) 



1=2 p=l+l 
X 



-^(p+l)...n-^i...p |(/ - 1) + ) (1 I ^(j,+iy..n^l...p 



(11 ^l-n^l-p b+) {l-\]^l-ni^l-p\{p + 1)" 

x(l-|^,..,^(,+i)...„|l+>' 



Similarly, from eq. (7.10) we obtain for the [g — s] contribution, 



Af\^a-X,...,n^)-i^^^^§0^y (7.39) 



where 



{{1-1)1) {l-\tl...nh\l^f 



fi"^ = -4 \^ M'' ~ v \-^ I -^^■■n^0 1-^ / 



K2 



=3 
X 



(7.40) 



(1 I ^(p+l)...„^0,i...p |(/ - 1) + ) (1 I i^{p+l)-n^4>,l-p 

{P (P+1)) 

{l~\^2-(l-l)^^,l-p {l-\^2-(l-l)^<p,l-p\{p + 1) + ) 

x(l-|^(,^,)...„^,,,..,|l+>^-^. (7.41) 

S<l>,l-p 

As indicated by eq. (2.18), and as explained in section VI, the amplitude with a fermion in the 
loop, An\(f), 1~, 2+, 3+, . . . , 71+), is equal and opposite in sign to the scalar-loop contribution 
A^n\(f), 1", 2+, 3+, . . . , n+) given in eq. (7.36). The soft Higgs limit of eqs. (7.36) and (7.39) 
onto the corresponding purc-QCD amplitudes (4.12) can easily be verified to obey eq. (5.13), 
just as in the quark case. 

Equation (7.39) completes the construction of all the finite one-loop QCD amplitudes 
containing a (p field plus multiple partons. 



VIII. RESULTS FOR TWO, THREE AND FOUR PARTONS 

In this section we collect explicit analytical results for up to four partons, using formulae 
from the previous section, as well as some results for divergent helicity configurations which 
have appeared elsewhere, in particular in refs. [31, 32]. 
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A. Complete results for two and three partons 



In appendix A, cqs. (A6), (A7), and (A8), we record the complete set of (pgg amplitudes, 
4'9^9^7 'Pg^g^: si'iid (j)g^g~ . The (j)qq amplitudes vanish (for massless quarks) by angular 
momentum conservation, so, rather trivially, we have the complete set of amplitudes for a 
single (f) plus 2 partons. 

The purpose of this subsection is to give the complete set of amphtudes for a single 
plus 3 partons. We start with the (t)ggg amplitudes. In the color decomposition (2.12), 
there is only one color structure, corresponding to the partial amplitude A^.i. First we 
record the finite helicity configurations, using our all-n results. Then, using these formulae 
and eq. (2.6), we decompose the one-loop i^gigfgi-amplitudes computed in ref. [31] into (pggg 
and (f>^ggg amplitudes. Finally we use parity to recast the results just in terms of (f)ggg 
amplitudes. These amplitudes will play a role in the recursive construction of one-loop 
Higgs amplitudes with more external partons. 

We express the one-loop results in terms of the corresponding Hggg or (l)^ggg tree-level 
amplitudes, 

4°^(i^,l^2^3+) = Af(0^l^2^3+) = ^py||^, (s.i) 

4°^(i7,l-,2+,3+) = Af(<^t,l-,2^3+) = -^l^lf^- (8.2) 

These amplitudes can be read off from eqs. (4.20), (4.22), and (4.28), with the help of 
eqs. (2.6) and (2.8). The [s] and [g — s] components of the finite loop amplitudes are, using 
eqs. (6.1), (6.2), (7.36), and (7.39), 

41(0,l+,2+,3+) = 0, (8.3) 

ASf-^l(0,1^2+,3+) = -2Af (0t,i+,2+,3+), (8.4) 
4^I(0,l-,2+,3+) = i^4°^(0t,i-,2+,3+), (8.5) 

4^-^1(0, l-,2+, 3+) = -2Af\<p\r,2\3+). (8.6) 
Using eq. (2.18), the full (f)ggg one-loop amplitudes are, 

>l3;i(0,l+,2+,3+) = -2 4'^)(0t,i+,2+,3+), (8.7) 



^3;i(</',l-,2+,3+) = 4°H</.^l-,2+,3+) 
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'^^3\' nJnJ si. 



(8.8) 



The one-loop Hggg amplitudes are [31], in our notation, 



As-AH, 1+, 2+, 3+) = Af\<l>\ 1+, 2+, 3+) L + ^ (l - 



Uf ns_\ £l2£23 + 523-531 + S31S12 

Nr NJ 



2\2 



A3;i(i/,l-,2+,3+) = 4°)(0t,l-,2+,3+) 
where 



2 

23 



(8.9) 
(8.10) 



C/3 = - 



In 



-S3I 



+ 



TT 



).(: 



-2 Li. 



2Li.(l-^U2Li2(l-^'l 



-S12 



+ 



-S23 



+ 



-S31 



£EUln(^^ln 

S3l\ 
kl) 



J (-S12 -S23\ J f-S23 -S3l\ J /-S3I -Sl2\ , , s 

with A;^ = S123, and fi is the scale originating from dimensional regularization. The one- mass 
box function Ls_i(ri,r2) is defined as [87], 



Ls_i(ri,r2) = Li2(l - ri) + Li2(l - r2) + lnrilnr2 - — . 



(8.12) 



We have extended the results of ref. [31] to nonzero values of Us by observing that the Uf 
and Ug terms are proportional to an amplitude for an off-shell gluon to split, via a fermion or 
scalar loop, into two on-shell, identical-hehcity gluons. (The opposite-hehcity case vanishes.) 
Such contributions are opposite in sign for fermions and scalars, as can be seen by contracting 
the off-shell gluon with an external qq pair, and then using a supersymmetry identity [76]. 

Using eq. (2.6), and subtracting eq. (8.7) from eq. (8.9), and eq. (8.8) from eq. (8.10), we 
have. 



A3;i(0t^l+,2+,3+) = A 



(O)r^t 1+ 9+ 



,^,1^2^3+) 
^3 + 2 + 1(1- 



?V , ^\ ■S12-S23 + -523-531 + -S3I-S12 



A3;l(0^1-,2+,3+) = Af\<p^,rXX)[U, + 2]. 



, (8.13) 
(8.14) 



The parity conjugates of these results constitute the remaining helicity amplitudes for plus 
three gluons, 

^3;i(0,l-,2-,3+) = 4°H0,l-,2-,3+)[[/3 + 2], (8.15) 
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^3;i(0,l-,2-,3-) = (0,l-,2-,3 



X 



nf Us \ S12S23 + S23S31 + S31S12 



2^2 



. (8.16) 



It is trivial to obtain the corresponding amplitudes for a pseudoscalar A plus three gluons 
via eq. (2.7), 



^3;i(Al+,2+,3+) = z4°)(0t,l+,2+,3+^ 

X 



A3;i(A,l-,2+,3+) = i4°)(0t,l-,2+,3+) 



nf ns \ S31S12 



^ 3V Nr NJ 



'23 J 



, (8.17) 
(8.18) 



Let us now turn to the (pqqg amplitudes. There is only one nonvanishing color structure 
in eq. (2.21), with partial amplitude ^3;!. We again express the one-loop results in terms of 
the corresponding Hqqg, or equivalently, (p^qqg tree amplitude. 



4°)(i/, 17, 2+ 3+) = 4^^(0^,17,2+3+) = - 



.[2 3]' 



(8.19) 



This amplitude can be obtained from eqs. (4.21) and (4.23), by using eqs. (2.6) and (2.8). 
We find from eq. (7.10), 

At'(^, ljX,3-)=2M + IMM = _24«.(^., 1- 2;, 3-) + iMM . (8.20) 
3 vs^> / > / > ; 2] 2 (2 3) ^ ' •'^ ' ' ^ 2 (2 3) ^ ^ 

The corresponding scalar amplitude ^3(0; l7,2j|r,3+) vanishes, but the R type amplitude 
does not, A§{(f), IJ, 2+, 3+) = -Af (0, 17, 3+, 2+) = -{i/2) (1 2) [2 3] / (2 3). Equation (2.25) 
then gives the full 4>qqg one-loop amplitude as. 



^3;i(0,l7,2+,3+) = -2Af(0t, 17,2+3+) + -(^1 + — 
= 4°)(0t,l7,2+ 3+ ' 



(1 2) [2 3] 
(2 3) 



(8.21) 



Prom the one-loop Hqqg amplitude [31] and eq. (8.21) we can extract the corresponding 
0^-contribution as above. The Hqqg amplitude is 



^3;i(//,17,2+,3+)=4°)(0t,l7,2+,3+) 



Vi + —:,V2 + -^Vs + -^V4 



7V2 



(8.22) 



with 



Vi = - 



V— S23^ V— Sai/J 6eV— Si9/ 



-S31 



13 / ji' 



6e V-S12. 
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83 _ 5k _ lsi2 
18 6 2 S23 



Vi = 



1 3i 



f fJ^'^ Y . J f-S23 -S3l\ . 7 Sr 1Si2 



_ 2 f n 



V. = - 



3eV- 



— S19/ 



1^4 = - - 



3e V— s 



12 



10 

8 

9 ■ 



(8.23) 

(8.24) 

(8.25) 
(8.26) 



Here Sr is a regularization-scheme dependent parameter, which fixes the number of hehcity 
states of the gluons running in the loop to (4 — 26 Re). For the 't Hooft-Veltman scheme [88] 
Sr = 1, while in the four-dimensional helicity (FDH) scheme [89, 90] Sr = 0. (Note that in 
the FDH scheme, the combination — V2 + V3 is equal to the function 11-^ appearing in the 
Hggg amplitude.) The Ug term can be deduced from the rif term using the fact that both 
arise from vacuum-polarization contributions. For this helicity configuration, the fermion 
and scalar loop contributions are not equal and opposite in sign. 
From eqs. (2.6), (8.21), and (8.22) then follows. 



A3;l(0^l7,2;,3+) = 4' 



(0), 



it,i-,2;,3^ 



2^J+A^^^+iV^^^ 



(8.27) 



The L, /?, / and s pieces of A3:i(0, , 2^ , 3~) can easily be extracted from the various color 
components of eq. (8.27) after applying parity, eq. (2.8). Note that the 512/(2523) terms in 
eq. (8.27) cancel similar terms in Vi and V2 in eqs. (8.23) and (8.24). 

Using eq. (2.7), the corresponding amplitude with a pseudoscalar A is found to be, 

1 



A3;i(A, 17, 2+ 3+) = I Af\4>\ 17, 2+ 3+) 



V^i + ^ + 4 + 
S23 



7V2 



V2 + 



S12 
S23 



Uf Us 



(8.28) 



B. Partial results for four partons 



The finite helicity amplitudes for a plus four gluons are given by. 



4^1(0,l+,2+,3+,4+) = 0, 



4^-^J(0,l+,2+,3+,4+) = -2Af (0t,i+,2 



aW(0,1-,2+,3+,4+) 



,3+, 4- 

v2 



;i2) [2 3] (13)" (14) [4 3] (13)^ 



(2 3)' (3 4) (41) (12) (2 3) (3 4)^ 



(8.29) 
(8.30) 
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(2 + 3) |4-) (3 + 4) |2-) 



Aj-^](0,l-,2+,3+,4+) = -2i 



(2 3) (3 4) S234 

(1-| (2 + 3) |4-)^ 



(8.31) 



(12) (2 3) (3-|(l + 2) |4-) S123 



(1- 


(3 + 4) 


2- 




' (3 4) (41) (3- 


(1 + 4) 


2~) Si34 



+ 



(^i)M2 4r 



[41] [12] (3-|(l + 2)|4-)(3-|(l + 4)|2-) s,^^ 
We have also computed the terms containing cuts in the following all-minus amplitude, 



32) 



At'\^, 1-, 2-, 3-, 4-) = (0, 1-, 2-, 3-, 4-)[C/4 + 2] , 



(8.33) 



where 



U, = — 



Ls_i 



/i 



'12- 
-S12 



+ 



"S123 —Si23' 



S34' 
—S23 



-S41 



-S34 



- Ls^i ( — —, — — ) - Ls_i 



"^234 —5234 ^ 
— S41 — S12 \ 



■S34I —5341/ ^—^412 —^412^ 

— Ls^°]"^(Si23, "S234; -523; ^^) " Ls^''j"^(s234, 5341; S34, /c^) 

-Ls2r(s341 kl) - Lsir (S412 

The two-mass box function Ls^™'^(s, t; mj, m|) is defined as [87], 



(8.34) 



2.e.„_2^2^ _ Li^fl-^VLi^A-^VLiJl-^ULi^fl-^ 



+ Li2 



2 2 



St 



t 

(8.35) 



This result agrees with the all-n cut-constructible expression An found previously in 
ref . [32] , when one substitutes n — A. As in the case of the Hggg amphtudes, the "2" added to 
C/4 in eq. (8.33) for >lj"'^(0, 1", 2", 3", 4") cancels against a "-2" in >l^"''((/)^ 1", 2", 3", 4") 
(the image under parity of ''(0, 1+, 2+, 3+, 4+) from eqs. (6.1) and (6.2)), when the 
and 0^ amplitudes are added together to produce A^^ ^\h, 1^,2^,3^,4^), or equivalently 
4'^'^^ of ref. [32]. 

The non-cut-constructible contribution to this amplitude was computed in ref. [32]. The 
result is 



A^r{(/),l-,2-,3-,4-) = S (1,2, 3, 4) + 5 (2,3,4,1) 

+ S- (3,4,1,2) + ^- (4,1,2,3), 



(8.36) 
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where 



"(1,2,3,4) = 



^13(4-1(1 + 3)12-)' (34)' J34)(41) 
[12f[2 3fsi23 [12f [12] [2 3] 

1 S12S34 + S123S234 - S?2 



(8.37) 



2 [12] [2 3] [34] [41] 

and the four terms in eq. (8.36) are the images of S (1, 2, 3, 4) under cyclic permutations 
of the gluon legs. This result was computed as a gluon and fermion loop contribution, but it 
can be extended to incorporate a scalar loop contribution using factorization properties. We 
checked the coUinear and multiparticle factorization limits of ALl~^\(j),l~ ^2>~ ^A~) and 
^^{(p, l-,2-,3-,4^) onto the corresponding (pggg and (pgg amplitudes. 

The soft Higgs limit of the scalar loop term ^4^(0, 1~ , 2~, 3-, 4~) in eq. (8.36) can easily 
be seen to obey eq. (5.13) with a factor of n_ = 4. The target pure-QCD amplitude, from 
eq. (4.5), is 

Equivalently, the Higgs amplitude ^^{H^ 1~, 2~, 3~, 4~) obeys eq. (5.3) for Z = 1 and n = 4, 
i.e. withafactorofn + 2Z-2 = 4 [32]. The [g-s] term A|f-^'(0, 1", 2", 3", 4") in eq. (8.33) 
vanishes in the soft hmit, because the tree amplitude 744*^^(1", 2~, 3", 4") vanishes, so that 
eqs. (5.11) and (5.12) for the soft limits of the divergent parts hold rather trivially. 

Finally, the amplitude for a pseudoscalar A plus four negative-helicity gluons is given via 
eq. (2.7) as 

^4;i(Al",2-,3-,4-) = -^A4;l(i/,l-,2-,3-,4-) + 2^A4(0^1-,2-,3-,4-) (8.39) 



—I 



<^(0,l-,2-,3-,4-)[[/4 + 4] 



(8.40) 



IX. CONCLUSIONS 



The production via gluon fusion of a Higgs boson in association with multiple jets provides 
an important background to Higgs production via electroweak vector boson fusion. Analytic 
representations of one-loop amplitudes containing a Higgs boson and four or more partons, 
interacting via the operator HiiG^vG^^, may be useful in quantifying this background. 
In this paper we have shown that on-shell recursive methods can be applied at the loop 
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level, to amplitudes containing a complex field plus multiple partons, where (f) = H + iA, 
with H a scalar Higgs boson and A a pseudoscalar. Previous applications of such loop-level 
relations were to pure QCD. There are three infinite sequences of finite one- loop amplitudes 
containing a single (p field and multiple QCD partons. We constructed recursion relations 
for all three sequences and provided the solutions. 

In particular, we presented compact formulae for the finite one-loop amplitudes with a 
single complex scalar 0, one quark pair, and n — 2 gluons of positive helicity, which are built 
from the primitive amplitudes y4^~*(0; j/) and ^4^(0; j/). Eqs. (7.10) and (7.31) provide the 
all-n forms of these primitive amplitudes with positive-helicity gluons. The corresponding 
tree-level quark-gluon amplitudes vanish, and hence these amplitudes are both infrared- 
and ultraviolet-finite. The amplitudes for a field and n gluons, all or all but one having 
positive hehcity, were presented respectively in eqs. (6.2) and (7.36). The case of one negative 
helicity was obtained by factorization from the finite quark-containing series of amplitudes. 
Together these amplitudes constitute all of the finite loop amplitudes for a single (f) field and 
multiple partons. Loop amplitudes containing additional quark pairs are always infrared 
and ultraviolet divergent, because the corresponding tree amplitudes are nonzero. 

Even though the corresponding tree-level amplitudes vanish, the finite one-loop ampli- 
tudes we computed enter next-to-leading order cross sections for Higgs-plus-jet production 
at hadron colliders, because the Higgs amplitude is a sum of (p and (p^ amplitudes. However, 
first the corresponding (f)'^ amplitudes have to be computed. Or equivalently (by parity), the 
(j) amplitudes for configurations with multiple negative-helicity gluons are required. 

We collected the complete set of one-loop (p amplitudes for two and three partons, as well 
as giving partial analytic results for four gluons, using also results from refs. [31, 32]. The 
four-quark Higgs results were presented analytically in cross-section form in ref. [21]. The (p 
amphtudes can be used to obtain amplitudes for a pseudoscalar A, as well as a scalar Higgs. 
We discussed the analytic properties of generic (p and H amplitudes in the soft Higgs limit, 
as the Higgs momentum vanishes. For the finite (p amplitudes this behavior is quite simple. 
The corresponding amplitude in pure QCD is recovered, multiplied by a factor of the number 
of negative helicities, n_. However, in the generic case the soft behavior is complicated, for 
many components of the amplitudes, by infrared divergences. 

The finite (p amplitudes also appear in factorization limits of the remaining, divergent one- 
loop (p amplitudes. Because of this property, their values will serve as input into a recursive 
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analytic construction of the latter amplitudes. We anticipate no conceptual problems in 
implementing such a unitarity-factorization bootstrap program [54, 56, 57], order-by-order 
in the number of negative-helicity gluons. 
Note added 

Soon after this work appeared, the first computation at NLO of Higgs production via 
gluon fusion in association with two jets was completed [91]. 
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APPENDIX A: NORMALIZATION OF y42;i(?5', 1+, 2+) 

In this appendix we give the result for the finite one-loop amplitude A2-^i{4>, 1"*", 2+) pro- 
duced by the effective Lagrangian (2.3), including the normalization factor, which feeds into 
all the other finite amplitudes discussed in the paper. We show that this normalization is 
consistent with results in the literature for the difference between the NLO QGD corrections 
to cross sections for inclusive production of a pseudoscalar Higgs boson, versus a scalar one, 
at hadron colliders. 

We have computed the one-loop helicity amplitudes for a scalar field H and pseudoscalar 
A, along with two positive-helicity gluons, A2;i{H,l~^ ,2~^) and A2;i(^, 1"*", 2+), whose nor- 
malization is defined according to eq. (2.12). In dimensional regularization, bubbles on 
external lines vanish. So there are no contributions from massless quark or squark loops 
because there is no direct coupling between the 4> field and massless fermions; in other words, 
^42;! has no dependence on Uf or n^. Accordingly (see eq. (2.16)), 




(Al) 




(A2) 
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There are only two nonvanishing diagrams in the scalar case: a triangle diagram, and a 
bubble diagram involving the four-gluon vertex in pure QCD. In the pseudoscalar case, the 
bubble diagram vanishes. In quoting the results, we do not perform any coupling renormal- 
ization, or renormalization of the effective Lagrangian operators in eq. (2.2). The results 
are. 



\—Si2/ L e 

2 



A;i(Al+,2+) = (-^) [-- + 4 + 0(6)1 (A 1^2+). 



(A3) 
(A4) 



The normalization of the leading poles in e agrees with the general structure of one-loop 

infrared divergences [86]. 

Because Af\(f),l^ ,2^) vanishes, and using eq. (4.20), we have A2°''(A, 1+, 2+) = 
i4°^(i/, l+,2+) = iAf\(j)\l+,2+). We can combine the H and A amplitudes (A3) and 
(A4) into the finite amplitude containing (f) — ^{H + iA), 

1 
2 

1/ 



A;i(0,l+,2+) = - A;i(i7,l+,2+)+zA.i(Al+,2+) 



+ 0{e) Af\H,l+,2+) 



+i 



-^ + 4 + 0(e) (//,!+, 2+) 



= -2Af{H,l+,2+)+0{e) 
= -24°^(0t,l+,2+) + O(e). 



(A5) 



Using eq. (2.6), we can obtain A2;i{4>\ 1^, 2+) by subtracting eq. (A5) from eq. (A3). We can 
then use parity, eq. (2.8), to get A2-i{(j),l'' ,2'). Now A2;i{(f),l^,2'^) vanishes by angular- 
momentum conservation. So the complete set of (pgg amplitudes is given, through 0{e), 

by 



A2;i(0,l+,2+) = -24°)(0t,i+,2+), 

A2;i(0,l^,2±) = 0, 



2 

A;i(0,i-,2-) = (j^y\-:^+2\Af\<p,i-,2-). 



2 

■^ + 2 



(A6) 

(A7) 
(A8) 



The finite one-loop amplitude (A6) also enters the difference in the NLO corrections to 
the inclusive hadron-collider production of a pseudoscalar vs. a scalar Higgs boson. In the 
large-mj limit, the part of the NLO corrections coming from real emission is identical in 
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the pseudoscalar and scalar cases [94, 95]. In terms of and (f)^ amplitudes, this result 
follows from the fact that the individual tree-level helicity amplitudes entering the real- 
emission subprocesses, gg Hg, qg — > Hq, gg — > Ag, qg — > Aq, plus processes related 
by crossing, are either MHV or anti-MHV — but not both. For every possible parton 
hehcity configuration, either the amplitude Af\(f), . . .) vanishes according to eq. (4.20) 
or eq. (4.21), or else the parity conjugate 0'^ amplitude Af\(/)\ . . .) vanishes. (Only when 
four partons enter the amplitude, do there begin to be helicity configurations for which the 
(j) and 0^ amplitudes are simultaneously nonvanishing. Such configurations, for example 
g~ g~ g~^ g~^ , are in the "overlap" of the tree- level MHV and anti-MHV towers [23].) 

If either ^3(0, . . .) or As{(f)\ . . .) always vanishes, then according to eqs. (2.6) and (2.7), 
the squared amplitudes for A must equal those for H, 



Af\A,l^X,^ 



± 



}±M2 



4°) (^,1^,2^,3^ 



± o± Q±M2 



(A9) 



and similarly for the qg — > Hq {qg Aq) process. Thus the real-emission contributions 
(normalized to the respective Born-level cross sections) are identical. 

The finite one-loop amplitude (A5) gives rise to a difference in the virtual corrections. 



(;^CrA,tree 



^^""^^ E.„.,|4°^(0, 1'^s 2^^^) + >if (0n I'^s s'^^) 

hi,h2 

-[^;l(0,l''S2'^^)+A.l(0t,l'^l,2''^)' 



3a. 



-4 Re 



A,,,{<f>, 1+, 2+)Af(0t, 1+, 2+)* + A2;i(0t, 1-, 2-)A'^\4>, 1-, 2-) 



Af\<P,l'\2'^-)-Af\<p\l'^\2'-] 



(0), 



47r 4o)(0t, l+, 2+)4°)(0^ 1+, 2+)* + 4°)(0, 1-, 2-)aP{<P, 1-, 2-)* 
An 

TT 



(-4)(-2) 



(AlO) 



However, eq. (AlO) does not take into account the difference in QCD corrections to the 
matching coefficients for the operators tr Gfj,,yG^'^ and tvG^u *G^^ (which is why wc put tildes 

over rfa-^(^)'i-'°°P). The operator trG^^G^'^ has coefficient C x [1 + lla,/(47r)] at NLO; the 
operator tr Gf^i,*G^'^ is unrenormalized. (The induced coupling of the pseudoscalar to the 
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divergence of the axial current does not contribute to the cross section at NLO, as discussed 
in section IV C.) Including this correction, we find 

^^yl,l-loop ^ Jf,l-loop ^~A,l-loop ^^H,l-loop II 

= 2 - 

(^(j A, tree da^'^'"''' rfcT^'*^'=<= (icT^'t''^^ 4?! 

^f6-ii)^ = ^. (All) 

This result agrees with the NLO calculations in refs. [94, 95, 96], which were also reconfirmed 
in the course of computing the NNLO corrections [97] . 



APPENDIX B: VALIDITY OF TREE-LEVEL ON-SHELL RECURSION RELA- 
TIONS WITH (p 

As outlined in section III, the derivation of the on-shell recursion relations relies on the 
vanishing of the amplitude A{z) as z ^ oo. For the purc-QCD case this result was shown 
in refs. [37, 40] . Below we extend the argumentation to amplitudes with a field. 



1. Vanishing of A{z — > oo) for the [— ,+) case 

In ref. [37] it was argued, based on the analysis of standard Feynman diagrams, that 
the amplitude A{z) in the purc-QCD case vanishes as 2; ^ oo for the case that the shifted 
helicities are hk = —, hi = +, with the corresponding shifted spinors given in eq. (3.1). 
The momenta are shifted according to eq. (3.2). Let us recall this argument: Any Feynman 
diagram contributing to a pure-gluon amplitude A{z) is linear in the polarization vectors of 
external gluons. Because only gluons k and I are shifted, only their polarization vectors can 
depend on z. These polarization vectors are given in spinor-helicity notation by 

[Xfj] ' (r/A) ' ^""'^ 

where 77, fj are fixed but arbitrary reference spinors. We see that these two polarization 
vectors behave asl/zas^;— >oofor the helicity configurations — —, hi — +, contributing 
a factor of z~'^ to the amplitude. The propagators each contribute a factor of 1/z, whereas 
triplc-gluon vertices are proportional to z in the z ^ 00 limit. Quartic vertices have no 
momentum factors. Therefore the worst possible contribution to A{z) stems from diagrams 
with only cubic vertices. A graph with (r + 1) cubic vertices has r propagators. Taking into 

54 



account the polarization vectors of the external gluons, the complete behavior as z — > oo is 
~ z^^^ I z^' X 2;^^. Thus the amplitude vanishes as I/2; as ;2 — > 00 for the helicity configuration 
^fe = -, = + ■ 

Now we extend this argument to the case where we attach a complex scalar to the 
amplitude, again with the above helicity configuration. Vertices with 3 and 4 gluons have 
the same 2; — > 00 behavior with and without a attached, thus the above argument applies 
directly to those Feynman graphs where the z-dependent momentum flows through only 
such vertices. However, the effective Lagrangian (2.3) also contains a vertex V^^g, coupling 
a field to two gluons. This vertex could potentially spoil the above power counting. 

Consider the set of graphs where the large 2;-dependent momentum flows through r triple- 
gluon vertices, a V^gg vertex, and r gluon propagators. The V^^g vertex is given by 

V7;r(^i,^2) = ■ k^T]^'^- - A^r^r +^£''^'^^''^"^^l.l^2.. , (B2) 

where k\ and k2 are the momenta of the two gluon lines. Naively, one might suspect such 
graphs to give nonvanishing contributions as 2; — > 00, because of the k^k^-X^Tva. in V^pp, 
which seems to behave as after the shift (3.2). However, the z^ contribution from the 
term ~ k^k^ vanishes using current conservation, after summing over the subset of graphs 
for which a fixed set of gluons attaches to each gluon line emanating from V^gg. That 
is, the two momenta ki and k2 are proportional to each other in the largc-2; limit, with a 
difference of order because k2 = —{ki + kr!)), and k^ carries no 2;-dependence. Thus we can 
replace k2^ {z)k'i^ {z) by —k'^^{z)k'^'^{z) in the large- 2; limit. The sum of graphs with a fixed 
set of on-shell external gluons attached to the momentum line ki forms a current J^^{z). 
Because current conservation, ki- Ji — 0, holds also for complex momenta shifted according 
to eq. (3.2), this sum of graphs vanishes as 2; ^ 00. The ^^-contribution from the first 
term in eq. (B2) vanishes on a graph-by-graph basis using eq. (3.2), and the 2;^-contribution 
from the last term vanishes on a graph-by-graph basis as well, using the antisymmetry of 
the Levi-Civita symbol. Thus the vertex behaves again at worst as z^, and the counting 
proceeds analogously to the pure-QCD case considered above, with the difference that a sum 
over certain subsets of graphs is required. In summary, tree-level amplitudes A{z) with and 
without (f) vanish at least as l/2;as2;— >oofor the helicity configuration hk — —, hi — + . 
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2. Vanishing of A{z — > oo) for the [+,+) and [— ,— ) cases 

The cases where the shifted gluons are both positive or both negative are shghtly more 
compUcated and rcqnirc a more involved argumentation along the lines of ref. [40]. In these 
cases, certain subsets of Feynman diagrams need to be summed to show the vanishing of 
A{z) as 2; — > oo already in the pure-QCD case. The extension to the ^-amplitude case works 
as in section B 1, summing over additional subsets of graphs. 

Here we recall the argument of ref. [40] , without going into too much detail, but pointing 
out the necessary modifications for the 0-case. We begin with the [+, +) case, and consider 
amplitudes which, besides gluons g^, g^, have also at least two gluons of negative helicity 
g'^,g^. Amplitudes with fewer than two negative- helicity gluons can be obtained via soft- 
gluon limits of such amplitudes. The argument is inductive in the number of external gluons. 
We proceed in two steps. First we perform a generalized shift on the spinors Aj, Ajt, A/, 

Aj — > Aj + z\k + zXi , Ajk ^ Ajfc — zXi , A; ^ A; — zXi . (B3) 

This triple shift preserves the sum of the momenta involved. Performing the same counting 
as in section B 1, the amplitude vanishes now as 1/2;^, or, without summing over the subset 
attaching to a V^^^ vertex, as 1/ z if the large- 2; momentum fiow hits such a vertex. In the 
latter case, however, for the argumentation below, we need to sum again over this subset as 
above, to obtain 1/ z"^ behavior. As usual, we obtain the recursion relation (3.6). 

Now we perform a second shift, this time only involving the positive-helicity gluons gf , g^ , 

\i + yXj , Xj Xj - yXi . (B4) 

If, in a given term of eq. (3.6), i and j are both contained in one of the amplitudes Al or Ar, 
then the entire ^/-dependence of the term is contained in one of these amplitudes; the other 
amplitude and the propagator arc independent of y. But Al and A^i have fewer legs than the 
full amplitude, and thus by the induction hypothesis the amplitude vanishes as y ^ 00. If i 
is on one side of the pole and j is on the other side, in a term of eq. (3.6), then the propagator, 
and therefore the value of the shift z at the pole, will obtain y-dependence. In the large-y 
limit, z y. By the above counting, the full amplitude behaves as 1/z^ ~ l/y^j where 
we have counted powers of internal vertices, propagators, and the polarization vectors of 
gluons i, /c, Additionally, with the shift (B4) we have to take into account the polarization 
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vector of gluon j, which now gives a factor of y because it has positive hehcity. In total, 
the amphtude now vanishes as l/y for large y. The argument for the [— , — ) helicity-case is 
analogous. 

To obtain this behavior, sums over two different subsets of graphs are required. In the 
pure-QCD case [40] , only the sum over the graphs on the left and on the right of the recursion 
relation (3.6) is necessary. For the case an additional sum over the subsets forming the 
current attaching to V^gg is required, if the large shift-dependent momentum flows through 
such a vertex. 

The above argument works if there are two extra gluons of opposite helicity to the shifted 
ones. If the amplitude under consideration does not have such a helicity configuration, then 
we can argue as in ref. [40] that we can obtain the desired configuration by taking the soft- 
gluon limit from an auxiliary amplitude with two extra gluons, as long as the auxiliary gluons 
are not adjacent to the marked gluons i and j. Moreover, in ref. [23] it was shown that the 
configurations A{(f), ±, +,..., +) vanish completely at tree level, cq. (4.20), without relying 
on on-shell recursion relations, and thus such an auxiliary construction is only necessary for 
the [— , — ) case. 
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